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SUMMARY 


An 11-month study contract, NAS3- 14422, was initiated to verify and increase 
the capability of the digital computer program developed from a previous contract 
NAS3- 13219. The objective of the present contract was to (1) experimentally 
verify the computer program developed from Contract NAS3-13219, (2) investigate 
different solution methods, (3) incorporate additional rotor dynamics parameters 
and (4) verify the expanded computer program. 

The final computer program simulates many important dynamic properties of a 
real rotor-bearing system. Included in the program are the effects of: 

1. Rotor slope coordinates 

2. Rotor material mechanical hysteresis in transverse shear and bending 
and torsional shear mode using viscous and/or Coulomb friction 
hysteresis coefficients. 

3. Torsional flexibility of rotor. 

4. Rotor transverse effects due to torsional and axial loading. 

5. Bearing in-phase and out-of-phase, anisotropic stiffness and damping 
force and moment coefficients. 

6. Bearing mass. 

7. Bearing transverse mass moment of inertia. 

8. Mount in-phase, anisotropic stiffness and damping moment coefficients. 

The program was written to analyze the general transient responses of nonsyn- 
chronous and nonaxisymmetric type rotor motion. The analysis of steady-state, 
synchronous and axisymmetric rotor motion can also be computed with the 
"HYSSTA" subroutine as a starting rotor-dynamic configuration. 

Evaluation of the digital computer computational speed of several mathematic 
approaches and integration methods were made. The validity of the computer 
program and accuracy of results was substantiated. 

Graphical (CRT) output capability is incorporated in the program. Detailed 
user's instructions are included in this report. 
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INTRODUCTION 


To further broaden the simulation capability and to improve the computational 
speed of the comptuer pgoram developed during contract' NAS3-13219, an 11-month 
contract (NAS3- 14422) was initiated. Extensions to the contract were subse- 
quently made to incorporate additional useful rotor dynamics effects. 

The major tasks accomplished in the contract were: 

1. Comparison of the initial computer program with the experimental results 
from Mark-25 test data was made. While only qualitative agreements 
were observed between the analytical and the test results, the computer 
results were compared accurately with those from other independently 
written computer program. 

2. Study of various mathematical approaches and integration techniques 
leading to a solution method which provides an optimum combination of 
computational speed and the predicted accuracy. 

3. Incorporation of the following additional useful rotor dynamics effects: 

a. Generation and application of the rotor slope influence coeffi- 
cients due to transient force and moment loading. 

b. Rotor mechanical hysteresis effects resulting from transverse 
shear strain. Viscous and/or Coulomb friction hysteresis coeffi- 
cients may be applied. 

c. Rotor transverse performance due to torsional and axial loading. 
(This automatically leads to the torsionally flexible rotor model.) 

d. In-phase and out-of-phase anisotropic bearing stiffness and damp- 
ing force and moment coefficients. 

e. Bearing mass and transverse mass moment of inertia. . 

f. In-phase anisotropic mount stiffness and damping and force and 
moment coefficients. 

/■ 

4. Final computer program verification. 

In Task (1), experimental verification of the computer program using Mark-25 
test data was made. The Mark-25 pump rotor is a very rigid rotor that has very 
small deflections and the calibration of the instrumentation could not be 
checked at all points due to the limited number of spot faces on the rotor. 
Therefore only some qualitative agreements between the data and analysis were 
observed. The computer program results have been compared with those from an 
independently written steady- state computer program and found to be valid and 
accurate. 
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As a result of Task (2) , the basic mathematical formulation for the computer 
program has been completely rewritten to accommodate the increase in new rotor 
dynamic parameters. Consequently, the size of the program has been substan- 
tially enlarged. 

The computer program which was written and developed on the G.E. time share 
computing system was converted to an IBM 360/370 version. Although substan- 
tial efforts in checking out the major operation of the program have been made, 
verification of all details of the program was not practical within the time 
and cost constraints. 

This final report presents the results of the computer program developed from 
1 March 1971 through 15 April 1973. 

The computer program, developed during this contract contains most of the major 
useful parameters believed to be required for close simulation of high speed 
flexible rotors. The program should provide a valuable tool for the analysis, 
design, and simulation of rotor-dynamic behavior. 

Rocketdyne strongly recommends a formal experimental verification of the com- 
puter program to further substantiate its validity and to gain confidence in 
the program use. 
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THEORY - MATHEMATICAL FORMULATION 


The governing mathematical formulation upon which the computer program was 
based is presented. There are two parts to the mathematical formulation; one 
is for the HYSSTA subroutine which generates the starting rotor dynamic config- 
uration, and the other is for the computation of the rotor dynamic performance 
through integration techniques used in subroutine FUND. Each part of the formu- 
lation will be separately described; The coordinate systems used are described 
in Figs. 1 through 4. 

I. MATHEMATICAL FORMULATION FOR SUBROUTINE HYSSTA 

The mathematical formulation is based on a steady-state, axisymmetric rotor 
motion of a multiple-station rotor-bearing system having axisymmetric rotor and 
bearing geometry, elastic modulii, force and moment stiffness and damping 
characteristics. The basic assumptions in the mathematical formulation are the 
linear elasticity and small deflections in the elastic rotor. The formulation 
considers the following parameters: 

1. Rotor transverse flexibility in shear and in bending. 

2. Rotor masses including transverse and polar mass moments of inertia at 
each rotor station. 

3. Rotor mass eccentricities and their corresponding phase angles. 

4. Rotor mass inertial misalignment angles and their corresponding phase 
angles. 

5. Two or more support bearings 

6. Linear rotor elastic modulii and bearing and other rotor support stiff- 
ness and damping characteristics. 

7. Bearing masses and transverse mass moments of inertia. 

8. Bearing in-phase and out-of-phase stiffness and damping force and 
moment coefficients. 

9. Mount in-phase stiffness and damping force and moment characteristics. 
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XYZ IS THE SECONDARY 
FRAME OF REFERENCE 
TO WHICH THE ROTOR 
FOUNDATION IS ATTACHED 
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ONLY TRANSLATORY MOTION 
BETWEEN THE TWO FRAMES 
OF REFERENCE IS PERMISSIBLE 



Figure 1. Relation Between the Secondary and an Inertial 
Frame of Reference 
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Figure 2. Rotor Transverse Coordinate Designation 
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Figure 3. Rotor Slpe Coordinate Designation 
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PHYSICAL DISTANCES AS SHOWN. THEY SHOULD BE INTERPRETED AS THE DISPLACEMENT 
FROM THEIR CORRESPONDING ZERO-LOAD POSITIONS. 



The following Eqs. (1) through (19) are used in HYSSTA and are solved simul- 
taneously to obtain a starting rotor dynamic configuration, using predetermined 
d>. values. 

T i 


<j>. I . + m.e 2 + I n . B 2 + 

1 I pi x l Di i J 

Vi [( Y i * *y) cos (♦i * °i) - ( X i * 8x) sin (*i * a i)_ 

l Dl B i[ 9 yzi cos (*i * Y i) - 9 xzi sin (*1 * Y i)]‘ 

! pi *1 B i [ 9 yzi sin {*i * Y i) * 9 xzi cos (*i * Y i)] - 0 


( 1 ) 
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F xi "A j X i * 8x - e i pi sin (*i * °i) *(hf cos (*i * “ 

+ |K. X. + C.X. + K .Y. + C . Y.l 
11 11 pi i pi l 

* [ K HDi Y i (h - hi hi) * C HDi Y i (V C Fi “Fi) 

‘ |(V*oi) [»Bi 


'BiK +B BiK (V X Bi +Y Bi' P BiK i 


+ K BiK 


C BiK(V X Bi +Y Bi" P BiK 


i H BiK 


+ D BiKlV X BiK + Y BiK P BiK ) + E BiK 


*81 


X EBxxi X Bi + X EBxyi Y Bi + ^EBxxi X Bi + ^EBxyi Y Bi 


F . = m . {y. + g + e. 
yi i | i 6 y 


( 2 ) 


t +i cos (*i * a i) ■ pi) 2 sin (» i * a i) | 

+ |K. Y. + C.Y. - K .X. - C .X. 

L i i ii pi i pi ij 

- [ K HDi x i (*i - *Fi W i) + C HDi X i (*i - C Fi “i)] 

[ N BiK +B BiK (\/ ( Bi +Y Bi _P BiK)] + K BiK J j C BiK (\/ C Bi +Y Bi" P BiK ) 

+ °BiK ^BiK + Y BiK “ P B1 k) + E BiK | Y Bi 


H BiK 


4, fc v -IT v 4 p v _ p v 

EBxxi Bi tBxyi A Bi fcBxxi Bi tBxy Bi 


( 3 ) 


where 


K EBxxi = 1/2 t K Bxxi + K Byyi> 
K EBxyi = 1/2 ( K Bxyi + K Byxi } 


(3a) 

(3b) 
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W ■ tc B „i * c Byyl) 

C EBxyi ' ( C Bxyi * Syxi J 

Xzi Wi * 'pj \ yi ij 


r- 




L*i siI> ( 

*i * V * (< 

y 2 cos (c 

>i * *i ) ] 


where 


* V 6xzl * '♦* 4 »‘ * V Vi + V Vi 

' Vl ft ■ >,) . , M ft . v ^ 

* >».» v„ • V, * W i„„ . w » w 

«• 

Myzi ' ' Di Vi ' V Vi h 

* 8 i ('pi ■ W pi «» (\ * yj ♦ (*.) 2 sln ^. + y jj 
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- [w Vi (»i - Vi V) * W e x2i (l t - c m %i jj 
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(3c) 

(3d) 


(4) 


(5) 


K EB4>xxi = V2 (K B0xx . 
K EB<f>xyi = V2 (K^. 
C £B(j>xxi = 1/2 (C B<J)xxi 

C EB(J)xyi a 1/2 (C^. 


B(J)yyi ) 

^Bcpyxi ^ 

+ C ■) 

B(f>yyi 1 

+ C •) 

B^yxi ' 


(5a) 

(5b) 

(5c) 

(5d) 
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n r 

s 

i=l 


Z.-S 


F . C. . + M . 

XI 1J XZl 


b ij] * Si * (Snb - Si) 4- - v 0 


( 6 ) 


n r 

St 


F . C. . + M . b. . 
yi i) yzi ij 
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bNB 
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- Y. 
J 
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n 
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|_ xi Fij xzi MijJ . L Bxzj 
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n 


Y -Y 
bNB bl 


E [f . T D . . + M . T . .1 +■ Ul 'y — - 0_ . = 0 
L Y 1 Fl J yzi MijJ L Byzj 


i=l 


C9) 


£ [(<* - Z i) F xi - M xzi] - 0 

i=l 


( 10 ) 


n 


E F . = 0 

xi 

i=l 


( 11 ) 


n 

ZJ (7q - Z.\ F . - M .1 = 

isl IV V yx yzij 


( 12 ) 


£ - o 


i=l 


yx 


(13) 


^EBxxi X Bi + K EBxyi Y Bi + C EBxxi X Bi + ^EBxyi Y Bi ” m Bi *Mi " ^Mi *Mi 
C Mi *Mi = 0 


(14) 
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• • M 

^EBxxi ^Bi ^EByxi *Bi + ^EByyi ^Bi ^EByxi ^Bi ” m Bi ^Mi " *Sli ^Mi 


C,,. Y u - = 0 
Ml Mi 


( 15 ) 


K 0 + k 0 +C 0*+C 0 

EB^xxi Bxzi EB<f>xyi Byzi EB<f>xxi Bxzi EBcftxyi Byzi 


_ _ • 

Bi Bxyi ^(j)Mi ^Mxzi ^(J>Mi ^Mxzi ^ 


^EB<J:yyi ^Byzi ^EBcf) i ^Bxzi ^EB4> i ^Byzi + ^EB4> i zi 


where 


^Bi ^Byzi *\f>Mi ^Myzi ^Myzi ® 


■Stt - 1/2 (*Mxi * Vi 5 

C Mi 1/2 ^ C Mxi + C Myi^ 

K.,.. = 1/2 (K... . + K, .. .) 
<j>Mi <j)Mxi cfiMyi 

C.,.. = 1/2 (C... . + C..„ .) 
<(>Mi ^ <J)Mxi <j)Myi 


Y. X. - X.Y. 

li li 

Pi ' X 2 + Y 2 
1 1 

0 .0 . - 0 .0 
yzi xzi xzi yzi 

“ Mi = e 2 ' . ♦ e 2 . 

xzi yzi 


(16) 


(17) 

(17a) 

(17b) 

(17c) 

(17d) 

(18) 
(19) 
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II. MATHEMATICAL FORMULATION FOR SUBROUTINE FUND 

The subroutine FUND is a time derivative generating subroutine based on the input 
parameters supplied by integration subroutines RKADAM. Thus the purpose of the 
mathematical formulation for FUND is to compute the corresponding accelerations 
from the displacements and velocities furnished by RKADAM. 

The formulation used in FUND is the most efficient one selected among various 
techniques investigated. It is a straight solution method without using simul- 
taneous equations approach that proved to be a more time-consuming method. 

The mathematical formulation used in FUND with its nomenclature discribed in 
Fig. 5 , consists of Eqs. (2) through (5) and (14) through (33). 


X i + 1 * X i 


Y i + 1 - Y i 


= £.0 . 
1 XZl 

/ o'.£. 
- 1 1 1 


| S x,i+1 

£ 2 

+ 1 
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\ 1 1 

3E.I. j 
1 

2E.I. 

i l 

= £.6 . 
l yzi 

/a’.£. 

I 1 i + 

\ 

S' 

£? 

+ 1 

l A'. G. 
\ i i 

3E.I. / 

li/ 

y.i+1 

2E.I. 

l i 


6 . - 0 . = 
XZ,1+1 XZl 


£? 

1 S' , . , + tt4- M’ 


£. 


2E.I. x,i+l E.I. xz,i+l 

l l * li ’ 


£ 2 


£. 


6 . , - 0 . = 
yz,i+l yzi 


f pi , i iii 

2E.I. y,i+l E.I. yz,i+l 

li 7 * li 1 * 


F . = S' ' - S' . 

XI X,l+1 XI 


F yi " S y,i+1 " S yi 


M . « - M' . , + M' . + S' . „ JL 

XZl XZ,1+1 XZl X,l+1 1 


M . a - M' . , + M' . + s' . , Z. 
yzi yz,i+l yzi y,i+l i 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 
(27) 
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Figure 5. Rotor Dynamic Configuration for the New Solution Method 
Used in the Computer Program 
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where 



C TVi * * /32 ' D oi 4 - D Ii 4 > U TVi A i 


(29) 

C TCi = YY12 (D oi " D Ii - 1 W TCi 


(30) 
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(32) 


M . = (0 . , - 0 .) T(}>RQ. 

Tyzi xz,i+l xzi r x i 


T «i - S-i <+i.l - *i> * C TVi «i.i - *i> * C TCi -,T 


- V 


l+i.l • h 


(33) 


The mathematical formulation used for the rotor material mechanical hysteresis 
effects due to in-plane (in-phase effects) and/or rotational (out-of-phase 
effects) rotor strain, is presented in Eqs. (34) through (47). Both the vis- 
cous and Coulomb friction models representing the hysteresis effects are in- 
cluded as denoted by subscripts V and C, respectively. 
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Equations (20) through (47) are used to solve for X^, Y^, 0 xz i and 9y Z j_ which 
are required for subroutine RKDAM. The <j>^ are determined from Eq. (28) using 

the computed values of X’i, Yi, 0 X zi anc * ®yzi- While not mathematical ly exact, 
this procedure for computing <f>i, Xj_, Y^ , 0 X zi an< ^ ®yzi has been demonstrated 
to give accurate results with substantial time saving over a simultaneous solu- 
tion method which is a mathematically rigorous procedure. 
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COMPUTER PROGRAM STUDY 
AND ITS VERIFICATION 


I. VERIFICATION OF THE EXISTING PROGRAM 

The computer program developed under contract NAS3- 13219 was to be verified by 
simulating a Mark- 25 liquid hydrogen pump rotor dynamic test performance. 

The simulation of the balancing data for the Mark 25 pump, without the inducer, 
was made for four combinations of speed and unbalance conditions. A comparison 
of the results with the test data indicates approximate agreement as shown in 
Tables I and II and Figs. : 6 through 9 . 

The unbalance configurations, designated as I and II, are defined in Table III; 
the Mark-25 rotor contour, including its balancing and the deflection stations, 
is described in Fig. 10 . . 

The comparison of test data with computer results indicates that substantial 
deviations exist at certain rotor stations. The degree of correlation varies 
among the four cases. With the exception of the deflection at rotor station 7 
for the unbalance configuration I, Fig. 6 , 7 and 9 indicate reasonable 
agreement in rotor deflection amplitude, while Fig. 7 shows approximate corre- 
lation in phase angles between the test data and computer results. 

To establish the validity of the computer program, a separate run from a pre- 
viously written steady-state rotor dynamic response computer program based on 
a different approach was made. The rotor deflections computed with the steady- 
state program are within 0.2% and the phase angles within 0.06 degrees of those 
computed with the present program. Hie steady-state program has been verified 
with other experimental data and found to be accurate. Thus, the discrepancies 
between the Mark- 25 test data and the computer results could be due to the 
inaccuracies of the test data. 

In examining the test data, particularly in Figs. 6 and 7 t it was observed 
that the rotor deflection at station 7 was substantially greater than the 
deflection at adjacent stations 5 and 12. This would not be realistic since 
with the large diameter portion of the rotor being sufficiently rigid, the 
deflection at station 7 should not be much more than that at station 5 or 12 
as demonstrated in the computer results. There are several possible causes for 
the inaccuracy of the test data. In reviewing the data tape, extensive DC 
drifts of the parameters were observed which may result in amplifier overloading 
and give rise to invalid scaled data. In the test data acquisition, a calibra- 
tion spotface on the rotor was provided only at one Bently transducer location 
near station 17. Hence, it is possible that calibration of the Bently trans- 
ducers at other stations may deviate from the standard value, thus causing 
erroneous scaling. 

It was also observed on the reduced test data that considerable noise and 
occasional variations in wave form existed. Difficulties have been experienced 
in deriving accurate rotor deflection and phase angle data where the magnitude 
of noise and wave irregularities were substantial in comparison with that of 
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TABLE I - COMPARISON OF TEST DATA WITH COMPUTER RESULTS (IN SI UNITS) 
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Table II - COMPARISON OF TEST DATA WITH COMPUTER RESULTS 
(SECONDARY SYSTEM; INCH AND DEGREE) 
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TABLE III - UNBLANCE CONFIGURATION (SI AND SECONDARY SYSTEM) 
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Figure 10. Mark- 25 Rotor (Without Inducer) Balancing and Deflection Stations 



the useful signal. Thus, it may be tentatively concluded that the discrepan- 
cies in the simulation are basically due to the inaccuracies accumulated in 
the process of test data acquisition. 

Experimental test data from Mechanical Technology Incorporated (MTI) was made 
available by NASA for an additional verification of the computer program. 

The test rotor-bearing configuration was modeled for the computer input. The 
rotor and bearing dimensions and properties used are indicated in Fig. 11 . 

The bearing stiffness and damping coefficients used were interpolated from the 
values at the speeds furnished by NASA as represented in Fig. 12 . 

The proximity transducer and the unbalancing plane locations are depicted in 
Fig. 13. A total of seven sets of test data was simulated with the computer 
analysis. The test data and corresponding unbalance and speed combinations are 
compiled in Table IV. The computer results, according to the rotor configura- 
tions and operating speeds, are shown in Figs. 14 through 20. Since the com- 
puter simulation involves only steady-state operation, a maximum of 0.001 
second of real time is used for the purpose of minimizing computation time. 

The "startup" rotor displacement coordinates are also included as a part of 
the computer outputs. 

Comparison between the computer results and the test data indicates that there 
are some similarities existing in the displacement and phase angle distribution 
along the rotor length, but the numerical values are not in close agreement 
as depicted in Figs. 14 through 20. 

To establish some confidence in the computer program, a parallel computer simu- 
lation of three sets of rotor test data using an independently written steady- 
state program was made. Close agreements were achieved between the results 
from the present transient program and those from the steady-state program. 


In an effort to improve the correlation between the computer results and the 
test data, computer simulations using other than the nominal bearing stiffness 
and damping coefficients were made. This action was taken assuming that the 
actual bearing stiffness and damping coefficients could deviate from the nom- 
inal computed values furnished by MTI. For simplicity purposes, the same 
constant ratio is maintained between the stiffness and damping coefficients as 
between the nominal stiffness and damping coefficients furnished by MTI. The 
results from this additional computer simulation (shown in Fig. 14 •) reveal that 
for some stiffness and damping coefficients used, some improved correlation in 
the magnitude of the displacement vectors and phase angles between the computer 
and test data was achieved. 

The correlation between the MTI test data and the computer results could prob- 
ably be substantially improved by using better defined values for some of the 
system constants. in the analytical model. 
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Figure 11. Rotor-Bearing Dimensions and Properties 




Figure 12. Bearing Stiffness and Damping Coefficients as Functions of Journal Speed 
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Figure 13. Proximity Transducer and Unbalance Plane Location 



TABLE IV. ROTOR TEST DATA AND RELATED UNBALANCE AND SPEED COMBINATIONS 
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Figure 14. Comparison of the Computer and Test Data for 
In-Line Unbalance at 6290 rpm 
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The aforementioned discrepancies between the computer results and test data 
could arise from one or more of the several reasons listed below: 

1. Accuracy of the test data. 

2. Accuracy of the unbalance data. 

3. Accuracy of the bearing stiffness and damping coefficients. 

4. Accuracy of the rotor modeling. 

5. Validity of the computer program. 

The first three items are parts of the test data package. No assessment of 
their accuracies can be made without adequate information concerning these data. 
The work regarding Items 4 and 5 was performed at Rocketdyne; therefore, their 
validity and accuracies will be discussed. 

Item 4 relates to the modeling of the rotor-bearing configuration. The major 
portion of the rotor modeling is sufficiently straightforward so that a definite 
procedure can be followed and accurate results can be predicted. Only at the 
locations where the rotor diameter changes abruptly, such as that at both sides 
of the central disk and at the locations of the bolted-on end disks, is special 
attention in modeling the rotor stiffness required. 

For the rotor sections with abrupt change in diameter, the local rotor stiff- 
ness is modeled according to that of the dotted lines shown in Fig. 11 to 
account for the bending stress distribution at the edges of the central rotor. 
These dotted lines represent the average stiffness of the shaft in the area of 
enlarged diameter. At the end disks of the rotor, the rotor section stiffness 
used corresponds to that of a slightly larger diameter shaft (as shown in 
dotted lines) than the nominal shaft diameter shown in Fig. 11. With the end 
disk-shaft configuration having a small shoulder and largely relieved shaft-to- 
hub contact surfaces, this shaft stiffness simulation appears to be reasonable. 
Furthermore, the end disks being the major end masses of the rotor, the effects 
of this shaft section stiffness on the rotor dynamic performance is much less 
than the case of a shaft section at the mid-span of a rotor. 

With reference to Item 5, "Validity of the Computer Program," the following may 
be added. A manual verification of the previous computer results for a 5- 
station rotor was made under contract NAS3- 13219 in Report NASA CR-72740. The 
computer results were found to be in precise agreement with that from the man- 
ual calculation. Additional verifications of the transient computer program, 
by comparing the results with those from the independently written steady-state 
rotor response program, was made once under contract NAS3- 13219 (Report NASA 
CR-72740) and twice (March 1971 and NASA Monthly Progress Report No. 10) under 
the present contract. The steady-state computer program, used in the verifica- 
tions, was based on a matrix iteration approach and it was shown to give good 
correlation with other test data. 
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Further verifications of the transient computer results including various damp- 
ing parameters were described in the October 1971 monthly report under the 
current contract. 

In conclusion, it may be stated that the transient flexible rotor computer pro- 
gram is believed to be valid and will yield accurate results. 

II. STUDY OF VARIOUS SOLUTION METHODS AND INTEGRATION TECHNIQUES 


To investigate the relative computational speeds and accuracies, various inte- 
gration techniques and computation methods were studied and results obtained. 

The solution methods studies are: 

1. Rotating coordinates using simultaneous equations without rotor slope 
parameters 

2. Stationary coordinates using simultaneous equations without rotor slope 
parameters 

3. - Stationary coordinates using matrix inversion without rotor slope 

parameters 

4. Stationary coordinates using independent equations including the rotor 
slope parameters 

The integration techniques used in combination with some' of the above solution 
methods are: 

1. Adams-Moulton predictor and corrector integration technique with 
variable time step 

2. Bulirsch-Stoer integration technique 

3. Three- time- level solution technique 

4. Fixed-step Adams-Moulton integration technique 

5. Fixed-step 4th order Runge-Kutta integration technique 

6. Fixed-step simple Euler integration technique 

7. Several modified, and some tolerance controlled Euler integration 
techniques 

The results in computation time, inverse of computation speed are presented in 
Fig. 21- through .27. The principles involved in the aforementioned solution 
methods and integration techniques are described as follows. 
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Figure 23. Computation Time Comparison for Steady-State 
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Figure 24. Computation Time Comparison for Accelerated 
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The solution methods are: 

1. Rotating coordinates using simultaneous equations without slope param- 
eters. This was the exact solution method developed under contract 
NAS3-13219. The acceleration of rotor dynamic displacements are solved 
in rotating coordinates at local rotor whirl frequencies. The purpose 
of using the rotating coordinates was to minimize the absolute magnitude 
of the second and third time derivatives of rotor displacements so to 
promote the computation speeds. This solution method while yielding 
very accurate computation accuracy, turned out to be a relatively slow 
process due to the necessary conversion required between rotating and 
stationary coordinates for each of the simultaneous solution procedures. 

2. Stationary coordinates using simultaneous equations without slope 
parameters. By removing the conversion process between stationary and 
rotating coordinates from the above rotating coordinates solution method, 
the computation process for comparable accuracy was found to be faster 
than using a rotating coordinate method for acceleration runs and slower 
for steady-state runs. 

3. Stationary coordinate method with matrix inversion and without rotor 

slope. parameters. Based on the mathematical formulation (Eqs. (2) •• 

through (13)), a set of inverted matrix coefficients using constant 4> 
per time step was computed according to the relationships below. 



(X; Y. X. Y. , 
v l i i x f 


<J>, <f> 40 


Y. = f . (X. Y. X. 'Y. , <J>, <J>, <j>) 
i yx i i x i* . T ’ 


<j> = f , (X. Y. X. Y X. Y. , <f>, <j>) 

T , <j) *• ,1 ■ l l , 1 l* T ' 


Thus the acceleration generating process, used for the rotor dynamics 
solution, involves only straight algebraic computation. Consequently, 
the process was substantially expedited. 


4. Stationary coordinates without simultaneous solution or matrix inversion, 
including the rotor slope parameters. The underlying principle of this 
method is first to compute the elastic and damping forces and moments 
due to relative displacements of the rotor, bearing and mount, and then 
apply these elastic damping forces and moments to compute the accelera- 
tion of the related masses and mass moments of inertia. This method was 
incorporated into the final computer program as it possesses an optimum 
.combination of computation speed and additional slope parameters. The 
’/.slope parameters improve the accuracy of computation results as compared 
with those of chord parameters developed under contract NAS3- 13219 . The 
mathematical formulation for this method is shown as Eqs. (32) through 
(39). ;; -v . v ' ■- . • • .• ; / ' • . • . 




Several integration techniques were used in combination with the aforementioned 
solution methods. The basic principle of each of the integration techniques 
are : 

1. Adams -Moulton predictor and corrector integration technique. * This 
integration method makes use of the third differences and a variable 
increment approach. The incremental time step is controlled by the 
specified input accuracy tolerance between the predicted and corrected 
solutions. The 4th order Runge-Kutta solution is used to initiate the 
computation. The predictor and corrector involves an iterative approach 
and is suitable for the area of rotor dynamic analysis where the magni- 
tude of second- and third-time derivatives vary widely. 

[ 2 ] 

2. Bulirsch-Stoer integration technique. This technique is based on 
a rational extrapolation in solving ordinary differential equations. 
According to a published reference, this technique has been demon- 
strated to be faster than some of the classic techniques. In applica- 
tion to the rotor dynamic analysis, it failed to show any advantage in 
computational speed over the Adams-Moulton predictor corrector method. 

r 3i 

3. Three-time-level solution method. The three-time-level (TTL) solu- 
tion method is based on the Taylor expansion of a Laplacian type of 
differential equation with time derivatives and constant terms. The 
method was first suggested by Dufort and Frankel and subsequently gen- 
eralized by Rocketdyne to solve multidimensional parabolic, hyperbolic, 
and elliptic differential equations in an application to transient 
thermal and related analyses. 

The application of the three-time-level technique has been shown to 
lead to less accurate results for a rotor without considering polar 
mass moment of inertia. For rotor including the polar inertia the 
technique caused computation .instability. Therefore, this method was 
discarded as an integration technique for the computer program. 

4. Fixed-Step Adams-Moulton technique. This technique uses the basic 
formulation applied in the Adams-Moulton predictor-corrector technique 
except that a fixed-time step is employed. It also uses the 4th order 
Runge-Kutta method as an integration starter. This method was incor- 
porated into the program as an option of the integration subroutine 
package which consists of steps 

a. Adams-Moulton predictor-corrector technique 


^ Hildebrand, F. B. , "Introduction to Numerical Analysis," McGraw-Hill, New 
York, 1956, pp. 199-201, 236-237. 

[ 2 ] 

Bulirsch, R. and Stoer, J. "Numerical Treatment of Ordinary Differential 
Equations ," Numerische Mathematik 8, 1-13 (1966). 
f3l 

L J Dufort, E. C. and Frankel, S. P., "Stability Conditions in the Numerical 
Treatment of Parabolic Differential Equations," Math. Tables Aids Comput., 
p. 135/152, 7 (1953). 
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b. Fixed-step Adams-Moulton technique 

c. Fixed- step Runge-Kutta technique 

5. Fixed- step fourth order Runge-Kutta integration technique. This is the 
classic fourth order Runge-Kutta technique using a fixed input time 
step and is incorporated as a part of the program integration package. 

6. Fixed-step simple Euler integration technique. This technique was 
studied only to demonstrate the limitation of a very simple integration 
technique against a common argument: that the computation accuracy should 
be achieved by applying a minimum time step but still results in compu- 
tation time saving due to the simple structure of the Euler integration 
technique. The computer results using this integration technique has 
shown extremely slow computational speed, for comparable accuracy with 
the other techniques evaluated such as Runge-Kutta and Adams-Moulton. 

7. Several modified and some tolerance controlled Euler integration 
techniques. The study results demonstrated that the simple Euler tech- 
nique requires longest computation time and a modified Euler provides 
shortest computation time among all the versions of Euler techniques 
investigated. The modified Euler is approximately 22 times faster than 
the simple Euler technique. 

Results of the study in computation speed for various combinations of solution 
methods and integration techniques are shown in Fig. 21 through 27. 

The results of the study have led to the selection of the following method and 
technique which give the best combination of computation and accuracy. 

1. Stationary coordinate solution method using independent equations and 
including rptor slope parameters. 

2.. Adams-Moulton predictor-corrector integration technique package which 
also includes fixed-step 4th order Runge-Kutta and fixed step Adams- 
Moulton integration technique. 

The overall improvement in computer time of the new method over that of the 
old one for the same program capability and comparable accuracy is substantial. 
The new computer program solution speed, including the rotor slope parameter, 
is approximately 2.3 times faster than the old program. Thecomputation-to- 
real time relationships for various rotor operating conditions are shown in 
Fig. 21 through. 26. 
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III. INCLUSION OF ROTOR DYNAMICS PARAMETERS 


During the contract period several significant and useful rotor dynamics param- 
eters were incorporated in the computer program as follows: 

1. Torsional flexibility of rotor 

2. Bearing mass 

3. Rotor material hysteresis 

4. Rotor transverse motion effects due to axial and torsional loading 

5. Bearing in-phase and out-of-phase anisotropic stiffness and damping 
force and moment coefficients 

6. Bearing transverse mass moment of inertia 

7. Mount in-phase anisotropic stiffness and damping moment characteristics 

The mathematical formulations for these parameters are described in this report 
under the section entitled "Theory - Mathematical Formulation". Computer dem- 
onstration of the effects of these included parameters are described in the 
subsequent sections. 

Rotor Torsional Flexibility Verification 

The rotor torsional flexibility contribution to the rotor spin drive torque, 
as defined in part (f) of Eq. (28) 

A drive torque i = (♦ i+1 “ < l > i ) - JCp (♦ i ”$ i _ 1 ) (55) 

has been demonstrated with the computer program. The computer results were 
validated with those from the hand calculation. 

Figures 28 and 29 contain the intermediate inputs and computer results for a 
3-station torsionally flexible rotor with a 106 in. -lb drive torque applied 
only at rotor station 2. No drive torques were applied to stations 1 and 3, 
which are torsionally connected to station 2 through the elastic shaft sections. 

The rotor configuration and properties used in the torsional flexibility com- 
putation are described in Fig. 30. 

Figures 31 and 32 present the dynamic performance of the torsionally flexible 
rotor. As shown in the figures, some torsional oscillation exists at all three 
stations. This torsional oscillation is induced from the suddenly applied 
drive torque in a nontorsional damping environment. 
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Figure 29. Rotor Drive Torque Step-Computation Results From 
Data Shown in Fig. 27 for a 106 In. -Lb Drive 
Torque Acting at Station 2 Only 
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Figure 31. A Computer Run Relative to the Data in Fig. 28 and 
Fig. 29 Using a 10^ In. -Lb Drive Torque Acting at 
Station 2 but With No Drive Torque at Stations 
1 and 3. 



THE AVERAGE REAL STEP-TIME FOR THIS PRINTOUT « 5.000E-0S SEC 

REAL TIME = 2.000E-03 SEC 

REVOLUTIONS ARRAY* 

8.026E+00 I.986E+00 2.026E+00 

SPIN SPEED ARRAY# RPM* 

7.974E+04 9.469E+04 7.974E*04 

ROTOR DISPLACEMENT VECTOR ARRAY# IN 

6* I 24E-03 2.769E-03 6.124E-03 

ROTOR VECTOR PHASE ANGLE ARRAY# DEGREES 
3.S21E+01 I.673E+02 3.521E+01 

ROTOR WHIRL/SPIN FREQ. RATIO ARRAY 

8»08IE-0l t.499E«-00 8.081E-0I 


THE AVERAGE REAL STEP-TIME FOR THIS PRINTOUT « 5.000E-05 SEC 

REAL TIME « 2.500E-03 SEC 

REVOLUTIONS ARRAY* 

2.732E+00 2.843E+00 2.732E+00 

SPIN SPEED ARRAY# RPM* 

1.046E+05 8.084E+04 1.046E+05 

ROTOR DISPLACEMENT VECTOR ARRAY# IN 

2.956E-03 1.897E-02 2.95«E-03' 

ROTOR VECTOR PHASE ANGLE ARRAY# DEGREES 
3.01IE+02 9.279E+OI 3.011E+02 

ROTOR WHIRL/SPIN FREQ. RATIO ARRAY 

1.127E+00 4.54SE-01 1.127E+00 


Figure 32. Continuation of Fig. 31 



Bearing Mass 


In verifying the GE computer program, a computation was made which included 
bearing mass, bearing and mount damping in a synchronous rotor motion. The 
magnitudes of the bearing mass and damping coefficients used in the computation 
were such that the computer results were distinct from those without these 
parameters thus allowing a comparison with accurate hand calculations for veri- 
fication. The related data used in the computations are, 


Bearing Mass = 

175.1 kg 

= 

1 (lb-sec^)/in. 

Rotor Spin Velocity = 

1000 radians/sec 



Bearing Stiffness = 

3.502 x 10 6 

Newtons/cm = 

2 x 

10^ lb/in. 

Bearing Damping 
Coefficient = 

1.751 x 10 4 

(Newton-sec) /cm = 

10 4 

(lb-sec)/in. 

Mount Stiffness = 

3.502 x 10 6 

Newtons /cm = 

2 x 

10^ lb/in. 

Mount Damping 
Coefficient 

1.751 x 10 4 

(Newton-sec) /cm = 

10 4 

(lb-sec)/in. 

Although the bearing damping coefficient is a major parameter 
the rotor dynamic configuration, it is, however, not directly 

in determining 
involved in the 


bearing mass equilibrium computation. Figure 33 gives the computer results 
from "STARTUP" computation. Verification of the configuration during the sub- 
sequent computation using an integration procedure was not made as it required 
a precise balanced torque to maintain a steady-state dynamic configuration. 
However, in a no-damping, steady-state run, the bearing mass coordinates com- 
puted from the subsequent integration procedure were found to be in agreement 
with those from the "STARTUP" computation as shown in Fig. 34. 

An axial projection of the bearing and mount dynamic configuration correspond- 
ing to the computed results presented in Fig. 33, is shown in Fig. 35. 

The individual force components and the force equilibrium conditions, as indi- 
cated by the values of relative error, are shown in Table V . 

The relative errors result predominantly from the computer printout roundoff 
errors and the equilibrium among the forces is thus established. 

From the results shown in Table V and Fig. 34, it may be concluded that the 
bearing mass parameter has been correctly formulated in the G.E. computer 
program. 
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BEARING MASS (LB-SEC 2 ) /IN 

O.OOOE+OO l.OOOE+OO O.OOOE+OO l.OOOE+OO O.OOOE+OO 

BCB, (LB-SSC)/IN ' • . 

0.000E-01 1 • 000E+ 04 0.000E-01 1.000E+04 0.000E-01 


THE COMPUTED STARTING ROTOR DEFLECTION COORDINATES ARE: 


X ARRAY, IN: 


-1..416E-03 

-1.3S8E-04 

-1.933E-03 

- 1 . 388 E- 04 

- 1 • 416E-03 

Y ARRAY, IN: 




1.68SE-04 

-3.^H6E-04 

3. 778 E- 04 

-3.4B6E-04 

1.685E-04 

VECTOR ARRAY, IN: 




1.426E-03 

3.752E-C4 

1.970E-03 

3.752E-04 

1.426E-03 

PHASE 

ANCLE ARRAY, 

DEGREES: 



I.732E+02 

2. ^63E+02 

1.689E+02 

2.4B3E+02 

1.732E+02 

XB(2>,Y«2>,XEC4>,YEC4>, IN 


• 


-6.724E-05 

-3. 549E-04 

-6.724E-05 

-3. 549E-04 



Note: BCB represents bearing damping coefficient. Only the 

2nd and lvth bearing mass and damping coefficient are 
used, in the computation. XB(2), YB(2), XB(4) and. YB(4) 
denote the rotor displacements from the bearing center 
in X- and Y-direction at the 2nd and 4th station 
respectively. 


Figure 33. Bearing Mass STARTUP Dynamic Configuration 
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BEARING MASS, (LB-SEC 2 )/lN 




O.OOOE+OO l.OOOE-Ol 

O.OOOE+OO 

1 .OOOE-Ol 

O.OOOE+OO 

BCB, (LB-SEC) /IN 

O.OOOE-Ol O.OOOE+OO 

0 • OOOE-0 1 

O.OOOE+OO 

0. OOOE-Ol 


THE COMPUTED STARTING ROTOR DEFLECTION COORDINATES ARE: 


X ARRAY , IN: 
4.998E-04 -1.479E-03 

y ARRAY, IN: 

1 • 000E-40 1 • 000S-40 

VECTOR ARRAY, IN: 
4* 998E-04 1 • 479E-03 

PHASE ANGLE ARRAY, 
1.800E+02 1 • 800E+02 


- 1 * 97 6E-0.3 


1 .000E-40 

1 .976E-03 
DEGREES: 

1 .800E+02 


-1 .479E-03 
1 .000E-40 
1 .479E-03 
1 .800E+02 


-4.998E-04 
1 .OOOE-40 
4.998E-04 
1 .800E+02 


XB(2>,YB(2)/ X3< 4 ) , YBC 4 ). HI 

-3.321E-04 O.OOOE+OO -3.321E-04 O.OOOE+OO 

MOUNT VECTORS: (2) AND (4), IN. 

1.147E-03 1.147E-03 

MOUNT VECTOR PHASE ANGLES: (2>'AND (4), DEGREES 

1 • 800E+02 1. 800E+02 


THE AVERAGE REAL 

STEP-TIME 

FOR THIS PRINTOUT = 5.000E- 

■05 SEC 

REAL TIME = 5. 

000E-03 SEC 




REVOLUTIONS = 

3.0002+00 

SPEED = 

3. 600000000E+04 

RPM 

VECTORS, IN? 

• 




4.994E-04 

1 .480E-03 

1 • 977E-03 

1 . 480E-03 

4. 994E-04 

PHASE, DEGREES: 





1 • 800E+02 

1 .800E+02 

1 . 800E+02 

1 • 800E+02 

1.8002+02 

WHIRL/SPIN FREQ. 

RATIO: 




9. 999E-01 

1 • OOOE+OO 

1 • 000E+00 

1. OOOE+OO 

9 . 999E-0 1 

MOUNT VECTORS: (2) AND <4>, 

IN. 



1.147E-03 

1 . 147E-03 





MOUNT VECTOR PHASE ANGLES: (2) AND (4), DEGREES 
1 * 800E+02 1 • 800S+02 


Note: Mount vector and mount phase angle represent the polar coordinates of the 

displaced bearing center. 


Figure 34. Bearing Mass Dynamic Configuration From Both "STARTUP" 
and Subsequent Integration Type of Computation 
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Figure 35. Axial Projection of the Bearing-Mount Dynamic Configuration 
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TABLE V - FORCE EQUILIBRIUM AT BEARING, MASS 


Force From 

X-Component 

Y-Coraponent 

Newtons 

(lb) 

Newtons 

(lb) 

Bearing 

Stiffness 

598.2 

2xl0 6 ) (0.6724xl0 -4 ) 
= 134. 48 

■ 

(2xl0 6 ) (3 . 549xl0~ 4 ) 
709.8 * 

Bearing 

Mass 

318.3 

1000) 2 (0.7156xl0 -4 ) 
= 71.56 

■ 

(1000) 2 (0.063xl0 -4 ) 
= 6.3 l 

Mount 

Stiffness 

636.6 

C2xl0 6 ) (0.7156xl0 -4 ; 
= 143.12 

56.} 

(2xl0 6 ) (0.063xl0" 4 ) 
= 12.6 | 

Mount 

Damping 

280.2 

•0 4 ) CIO 3 ) (0.063x10 

= 63 . 

a 

3183 . | 

LO 4 ) CIO 3 ) (0 . 7156xl0 -4 
= 715.6 | 

Relative 

Error 

4 

^Force 

~1 J , , ~ 0.04% 

0 . 5 £ Force J 

1 

4 ’ 

J] Force 

7 J - 0.07% 

0,5 2 | Force| 
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Hysteresis 


In the verification of the hysteresis formulation in the computer program, a 
hysteresis moment computation was made by using an artificial value of 4 >-uj c ^. 
The resulting hysteresis moments generated has been verified to be correct. 

Rotor Transverse Motion Effects 


In both the axial and torsional loading computations, a 3-station symmetric 
rotor-bearing configuration was used for each loading. The two configurations, 
shown in Fig. 36 , were identical except for different polar mass moments of 
inertia. No basic reasons for using different rotor configurations exist, they 
simply represent alternate rotor designs. 

Axial Loading Effects on Rotor Transverse Motion . In the computer analysis, 
various combinations of axial loading, rotor spin speed, and damping parameter 
were included in the computer runs, as presented in Table VI . 

Three 1-rpm rotor spin speed runs (Runs 1 through 3) were first made to simulate 
the static performance of an axially loaded, hinged-end, uniform cross-section 
column having an initial elastic deflection. In Run 1, a 44,482 Newtons (10 4 
lb) axial compressive loading was used. The midspan (Rotor Station 2) rotor 
deflections are depicted in Fig. 37. In Run 2, the same operating conditions 
were maintained except that a 44,482 Newtons axial tensile loading was applied, 
instead of the compressive one. The corresponding computer results are shown 
in Fig. 38 , which indicates a reduction of rotor deflections leading to stable 
operation as expected. Run 3 is similar to Run 1, except that a smaller axial 
loading of 222 Newtons (50 pounds) was used. The results are presented in the 
lower portion of Fig. 38. Although the rotor deflection magnitude in Run 3 is 
much reduced from that in Run 1, the general trend of monotonically increasing 
deflection remains the same as in Run 1 as shown in Fig. 37. 

The critical buckling load for the hinged-end, uniform cross-section column 
was computed from the Euler formulation (P C ritical = ^ EI/L.2) to be 259 
Newtons (58 pounds). The reasons for the instability trend observed in Run 3 
with below the critical compressive loading are twofold; first, a sectional 
linear moment axial loading model was used in the analysis, as compared with 
the actual nonlinear moment-axial loading function, and second, the inherent 
instability effects of the mass-unbalance rotor dynamic loading considered in 
the analysis. For a rotor having a reasonable number of stations, the computer 
analysis will closely approximate the exact non-linear moment-axial loading 
function of an axially loaded rotor. 

In Run 4 through 9, a moderately high rotor spin speed of 36,000 rpm in combin- 
ation with different axial loadings and damping coefficients was used. The 
rotor deflection results for these runs are depicted in Fig. 39 through 41. 

A comparison study of the undamped axial loading results (Runs 1 through 5 and 
run 8) suggests that a high rotor spin speed tends to resist buckling instabil- 
ity even under large axial compressive loading. While a larger than critical 
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ROTOR STATION NO. 


1 


2 


3 


MASS: Kg (LB) 

22.68 

(50) 

22.68 

(50) 

22.68 

(50) 

TILIN' S VERSE INERTIA 

Kg-Ctf^LB-IN 2 ) 

176 

(60) 

176 

(60) 

176 

(60) 

POLAR INERT LA 

Kg- CM* (LB- IN* ) 

FOR AXIAL LOADING 

146 

(50) 

146 

(50) . 

146 

(50) 

FOR TORSIONAL LOADING 205 

(70) 

205 

(70) 

205 

(70) 


MASS ECCENTRICITY CM (IN).[ 0.000762 0.00254 0.000762 

AT 45 DEG PHASE ANGLE (0.0003) (0.001) (0.0003) 

INERTIA MISALIGNMENT 0 0-0 



ROTOR MATERIAL PROPERTIES: 
DENSITY 

YOUNG MODULUS OF ELASTICITY 
SHEAR MODULUS OF RIGIDITY 
POISSON'S RATIO 


0 

2.068X10 7 NEWT0N/CM 2 (3 x 1G 7 PSI) 
0.793XI0 7 NBfTON'/CM 2 ( 1 . 1 D x 1 0 7 PS I) 

0.3 


Figure 36. Rotor Configuration 


66 



PARAMETER INCLUDED 


RUN 

NO. 

ROTOR 
SPIN SPEED 

AXIAL LOADING 

DISPLACEMENT 

DAMPING 

COEFFICIENT 


RIM 

NEWTONS 

(LB) 

NEWTONS -SBC/cm 
(lb-sec/in) 

i 

i 

1 ! 

1 

44,482 *C 

(10 4 c) 

0 

2 

1 

44,482 *T 
(10 4 T) 

0 

3 

1 

222. C 
(50 C) 

0 

4 

36,000 

222. C 
' (50 C) 

0 

5 

36,000 

44,482 C 
(10 4 c) 

0 

6 

36,000 

44,482 C 
(10 4 C) 

1,751. 

(10 3 ) 

7 

36,000 

444,822 C 
(lO 5 C) 

1,751 

(103) 

8 

36,000 

444,822 C 

(105 c) 

0 
























































loading would definitely result in rotor instability under a non- rotating, 
static condition, a high-speed rotor system of the same design may operate 
without exhibiting an unstable trend. The mechanism for the apparently more 
stable motion with a high speed rotor under compressive axial loading appears 
to be similar to that causing the stable motion for a super-resonant rotor 
speed operation. In such a super- resonant rotor speed operation, the mass 
acceleration stiffness exceeds that of the combined rotor-bearing stiffness 
which would lead to instability if not for the dynamic effects of rotation. 

Runs 6 and 7 included viscous damping in the rotor-bearing system under a large 
axial compressive load. The corresponding rotor deflections are plotted in 
Fig. 40, which indicates a substantial reduction in rotor deflection and oscil- 
lation compare with those in Run 5, Fig. 39, and Run 8 , Fig. 41, where damping 
was absent. '■ 

Figure 41 depicts the results from Run 8 which was conducted under the same 
operating conditions as Run 7, except without damping. In Fig. 41, Run 8 ex- 
hibits a trend of larger rotor deflection over that in Run 5 (Fig. 39) due to 
larger axial loading and a much stronger oscillation as compared with Run 7 due 
to the absence of damping. 

A review of the rotor dynamic deflection performance of the eight computer runs 
suggests that the rotor rotation, particularly at high rotational speed, and 
damping axe stabilizing factors against the buckling effects of an axial com- 
pressive loading. The stabilizing effects include the minimization of deflec- 
tions as well as oscillations induced by axial loading. 

Torsional Loading Effects on Rotor Transverse Motion. In the analytical demon - 
stration of the torsional loading effects, eight computer runs (Runs 9 through 
16), which include various combinations of rotor speed, torsional loading and 
damping coefficient, were made as stated in Table VII. The computer results 
are depicted in Fig. 42 through 44. The rotor bearing configuration used in 
this analysis is defined in Fig. 36. 

For each of the computer runs, positive drive torque was applied at Rotor 
Station 3 (rotor right end) balanced with an equal and opposite drive torque 
at Station 1 (rotor left end). Thus* the rotor spin acceleration at the mid- 
rotor station (Station 2) will remain zero and the average speed at the ends 
of rotor will be that of the initial starting rotor speed. Using the opposite 
drive torques at the ends of the rotor, the rotor transverse effects due to 
torsion without those due to the average rotor acceleration caused by unbalanced 
torque application may be clearly observed. 

To simulate the torsional loading effects under the nonrotating, static condi- 
tions, two 1-rpm and two 10-rpm runs were included; Runs 9 and 10 are those at 
1-rpm and Runs 11 and 12 were made at 10-rpm. The computed rotor deflections 
at rotor ends (Stations 1 and 3) or the rotor deflections at Station 1 when 
they approximate those at Station 3, are depicted in Fig. 42 through 44 , in the 
torsional loading analysis. 
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TABLE VII - TORSIONAL LOADING RUNS 


RUN 

NO. 

PARAMETERS INCLUDED 

ROTOR 
SPIN SPEED 

RPM 

♦BALANCED 
DRIVE TORQUE 

NEWTON-CM 

(LB-IN) 

DISPLACEMENT DAMP- 
ING COEFFICIENT 

NEWTON-SEC/ CM 

(lb-sec/in) 

1 

SLOPE DAMPING 
COEFFICIENT 

NEWTON-CM-SEC/RAD LAN 

(lb-in-sec/radian) 

9 

1 

7909 

(700) 

0 

0 

10 

1 

i 

11.2985x10° 

(10 6 ) 

0 

0 

i 

11 

1 

i 

10 

7909 

(700) 

0 

0 

12 

■ 10 

1 

11.2985x10° 

( 10 6 ) 

0 

0 

13 

36,000 

■ 

1 

7909 

(700) 

0 

0 

14 

1 

36,000 

j 

1 

11.2985x10° 

( 10 6 ) 

0 

0 

15 

! 

i 

36,000 

11.2985x10° 

do 6 ) 

17513 

( 10 4 ) 

0 

16 

36,000 

11. 2985x1 0 6 

do 6 ) 

0 

11.2985xl0 4 

( 10 4 ) 


^Positive Drive Torque at Station 3 and An Equal and Opposite Drive Torque at 
Station 1. 
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INCH . ROTOR DEFLECTION AT STATION 1, 10 INCH 















ROTOR DEFLATION AT STATIONS 1 AND 3, 10~* INCH 



Figure 43. Torsional Loading Effects at 36,000 rpm. 
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RUN 15 WITH 17,515 NEWTON-SEC/ CM DISPLACEMENT DAMPING COEFFICIENT 
AND 11,2985x10" NEWTON-CM DRIVE TORQUE. i — : — t — — r— r 


i iaBaBBSMB Siia ja l 

:ssnssassss:u: 1 


RUN 16 WITH 11 .2985x10 i NEWTON-CM-SEC/RADIAN SLOPE DAMP 
COEFFICIENT £ND 11.2985x10° NEWTON-CM DRIVE TORQUE, ~ 4 


IBS 


RUN 16, 11 •’2985x1 0° NEHTON-CM DRIVE TORQUE WITH 
11-2085x10^ NEWTON— CM— SEC /R ATTT AM ST CPE TUMPING CHEFF 


TIME, SECOND 

Figure 44. Torsional Loading Effects Including Damping at 36,000 rpm 








The midrotor deflections remain basically constant in all the runs under the 
balanced torque conditions. In Fig. 42 large oscillations are observed for the 
11.2985 x lOo Newton-cm (10^ lb-in.) torque runs, while for a 7909 Newton-cm 
(700 lb-in.) torque loading, the oscillation in rotor deflections are drastic- 
ally reduced. The effects of a small change of rotor speed from 1 to 10 rpm is 
minimal as expected. The results presented in the upper portion of Fig. 42 
suggest a trend of instability, while the lower portion may lead to stable 
operation. 

For the static, nonrotating, uniform cross-section rotor, the critical torsional 
buckling load is computed to be 8363 Newton-cm (740 lb-in.) from the classic 
formulation; torque cr ^ t ^ cal = 2 tt EI/L. The value of smaller torsional loading 
(7909 Newton-cm) used in the analysis is below the static critical buckling 
loading, where the larger torque loading far exceeds the critical value. 

Runs 13 through 16 were made at a moderately high operating speed of 36,000 rpm. 
A lower torsional loading of 7909 Newton-cm was used in Run 13 and the higher 
loading of 11.2985 x 10^ Newton-cm was applied in Run 14. No rotor damping was 
included in Runs 13 and 14. The computer results are depicted in Fig. 43 . Run 
13, shown in the upper portion of Fig. 43, indicates a high degree of stability, 
while Run 14, which used large torsional loading, shows oscillations and a di- 
vergence in rotor deflection versus time relationship. 

Runs 15 and 16 included displacement and slope damping, respectively. In both 
the runs the high torsional loading of 11.2985 x 10^ Newton-cm was used. Rotor 
deflection data for Run 15 are depicted in the lower portion of Fig. 43 which 
indicates a modification of the deflection versus time characteristic from a 
divergently oscillatory in the no-damping case of Run 14 to a monotonically 
increasing behavior, except near zero time point. The explanation for the per- 
formance in Run 15 versus Run 14 is that although the oscillation energy has 
been dissipated through damping, the divergent, unstable factor induced by the 
super-critical buckling torsion still persists. 

Run 16 was computed under the effects of slope damping which acts against the 
change of slopes, but not against the change of deflection. Clearly no notice- 
able damping effects on rotor deflections are observed, as shown in the lower 
portion of Fig.' 44. However, a comparison of the rotor slope versus time func- 
tions between Run 15 and 16, as shown in the upper portion of Fig. 44, indicates 
that the slope damping was demonstrated to be substantially effective in mini- 
mizing the slope oscillation in Run 16. 

The study of the torsional loading computer results revealed that the applica- 
tion of a below critical buckling torque would not cause rotor oscillation or 
instability in the high speed rotor operation, particularly in the presence of 
rotor damping. High speed rotation and damping parameters are, therefore, 
considered to have a substantial stabilizing effect for rotors under torsional 
loading. 
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Bearing and Mount Parameters 


These added parameters include: 

1. Bearing in-phase and out-of-phase stiffness and damping forces, 

2. Bearing in-phase and out-of-phase stiffness and damping moments, 

3. Bearing transverse mass moments of inertia, and 

4. Mount in-phase stiffness and damping moments. 

To demonstrate the validity of the GE coding of these parameters, sample rotor 
dynamic performance computations were made. Figures 45 through 53 represent 
the rotor dynamic performance according to the input parameters defined in 
Fig. 54. Each of the computer runs consists of two parts; part a is the startup 
rotor dynamic configuration generated in subroutine STARUP and part b repre- 
sents the integration results computed in subroutine RKADAM. The agreements 
between the rotor performance from STARUP and RKADAM of each of the computer 
runs indicate that the coding in RKADAM is consistent with that in STARUP. The 
validity of the startup rotor configuration has been demonstrated in three ways: 

1. Force and moment equilibrium among those due to the displacements 
and slopes of the bearings and mounts and bearing masses and mass 
moments of inertia as follows: 

Bearing force + bearing mass force = mount force (along X and Y axes) 

Bearing moment + bearing inertia moment = mount moment (in X-Z and 
Y-Z planes) 

2. Displacement and slope consistency among those of the rotor, bearings 
and mounts (i.e., the vectorial sum of displacements or slopes of 
bearings and mounts equal those of the corresponding rotor displace- 
ments or rotations). 

3. The bearing and mount forces and moments correspond to the products 
of their displacements and slopes, their first time derivatives and 
the stiffness and damping coefficients. 

Since there is a unique solution to a rotor system at a set of operating condi- 
tions, other than those unstable rotor dynamic configurations due to nonlinear 
bearing stiffness, etc., the consistent solutions represented in Figures 45 
through 53 are the only valid ones. 

In the aforementioned computer runs, axisymmetric bearing and mount stiffness 
and damping coefficients were used although the general nonisotropic coefficient 
capability is included in the RKADAM subroutine. The reason for the use of 
axisymmetric coefficients was to provide an exact comparison of rotor perform- 
ance between that generated from RKADAM and that from STARUP which is limited 
to axisymmetric coefficients as planned. 

The other rotor design and operating conditions which are not contained in Fig. 
36, but common to all the 9 computer runs, are described in Fig. 55. 
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THE COMPUTED STARTING KOTOR DEFLECTION COORDINATES ARE: 
VECTOR ARRAY* IN: 

1.1409E-04 9.3123E-05 1.1409E-04 

PHASE ANGLE ARRAY* DEGREES: 

2.2500E+02 2.2500E+02 2.2500E+02 

BEARING DISPLACEMENT VECTOR* IN 
1 • 3801 E- 04 1.380IE-04 

BEARING DISPLACEMENT PHASE ANGLE* DEGREES 
2.2500E+G2 2.2500E+02 ' 

MOUNT VECTORS* IN 

2.3918E-05 2.391SE-05 

MOUNT VECTOR PHASE ANGLES* DEGREES 
4.5000E+01 4. 5000E+0 I 

ROTOR SLOPE VECTORS 

1 •2827E-04 9.5349E-I4 1.2827E-04 

ROTOR SLOPE VECTOR PHASE ANGLE* DEGREES 
2.2500E+02 2.2060E+02 4.5000E+01 

BEARING SLOPE VECTORS 

O.OOOOE+OO O.OOOOE+OO 
BEARING SLOPE PHASE ANGLE * DEGREES 
4.5000E+0I 4 • 5000E+0 I 

MOUNT SLOPE VECTORS 

1.2827E-04- 1.2R27E-04 

MOUNT SLOPE PHASE ANGLE* DEGREES 
2* 250QE+02 4.5000E+01 . 

BEARINu ru.iioc. IN X- DIRECT I ON 

-9.7590E+01 -9 • 7 590E+0 1 

BEARING FORCE IN Y-DIRECTION 
-9.7590E+01 -9.7590E+01 

BEARING MOMENT IN X-Z PLANE 
O.OOOOE+OO O.OOOOE+OO 

BEARING MOMENT IN Y-Z PLANE 
O.OOOOE+OO O.OOOOE+OO 

M01WT FORCE IN X-DIRECTI0N 
3* 3825E + 0 1 3.3S25E+01 

MOUNT FORCE IN Y-DIRECTION 
3.3825E+01 3. 3825E+01 

MOUNT MOMENT IN X-Z PLANE 

O.OOOOE+OO O.OOOOE+OO 

MOUNT MOMENT IN Y-Z PLANE 

O.OOOOE+OO O.OOOOE+OO 

BEARING MASS FORCE IN X-DIRECTION 
1.3142E+02 1 .3 i 42E+02 

BEARING MASS FORCE IN Y-DIRECTION 
1.3I42E+02 1.31 42E+02 

BEARING INERTIA MOMENT IN X-Z PLANE 
O.OOOOE+OO O.OOOOE+OO 

BEARING INERTIA MOMENT IN Y-Z PLANE 
O.OOOOE+OO O.OOOOE+OO 


NOTE: FOR FIG. 45 THROUGH 53 

THE UNITS USED IN THE COMPUTER 
RUN ARE: 

DISPLACEMENT = INCHES 

SLOPE = RADIANS 

FORCE = LB 

MOMENT = IN- LB 


Figure 45. Computed Starting Rotor Deflection Coordinates 
Run 2 



THE AVERAGE REAL STEP-TIME FOR THIS PRINTOUT = 2.500E-05 SEC; 

REAL TIME = 1.000E-03 SEC 

REVOLUTIONS ARRAY: 

I • S9 1 SE+OO 1.5915E+00 1.5915E+00 

SPIN SPEED ARRAY# RPM : 

9» 5493E+04 9.5493E+04 9.5493E+04 

ROTOR DISPLACEMENT VECTOR ARRAY# IN 

I • 1 409E-04 9.3I24E-05 1.1409E-04 

ROTOR VECTOR PHASE ANGLE ARRAY# DEGREES 
7.7958E+0! 7.7959E+01 7.7958E+01 

BEARING DISPLACEMENT VECTORS 
1.380IE-04 1 * 380 1 E- 04. 

BEARING DISPLACEMENT PHASE ANGLES# DEGREES 
7«7958E+01 7.7958E+01 

MOUNT DISPLACEMENT VECTOR ARRAY# IN 
2.39I4E-05 2.391-aE-ns 

MOUNT VECTOR P!'- SE ANC ( E DEGREES 

2.5796E+02 2.5796E+02 

ROTOR WHIRL/SPIN FREQ. RATIO ARRAY 

1 • OOOOE+OO 1 •OOOOE+OO I.OOOOE+OO 

ROTOR SLOPE VECTORS 

1.2825E-04 1.1498E-13 1..2825E-04 

ROTOR SLOPE PHASE ANGLES# DEGREES 

7.7958E+01 3.4306E+02 2.5796E+02 

BEARING SLOPE VECTORS 

3» 3399E- 1 5 2.2888E-I5 

BEARING SLOPE PHASE ANGLES# DEGREES 
1 • 9S42E+02 3.37I7E+02 

MOUNT SLOPE VECTORS . 

1.2825E-04 1 • 282SE- 04 

MOUNT SLOPE VECTOR PHASE ANGLES# DEGREES 
7.7958E+0I 2.5796E+02 


Figure 45, Concluded 
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THE COMPUTED STARTING R0T0R DEFLECTION COORDINATES AREt 
VECTOR ARRAY* INl 

1.0167E-04 9.3438E-05 1.0167E-04 

PHASE ANGLE ARRAY* DEGREESl 
2.2072E+02 2.2512E+02 2.2072E+02 

BEARING DISPLACEMENT VECTOR* IN 
I.0013E-04 1.0018E-04 

BEARING DISPLACEMENT PHASE ANGLE* DEGREES 
2.1088E+02 2. 1088E+02. 

MOUNT VECTORS* IN 

I. 7362E-05 I.7362E-0S 

MOUNT VECTOR PHASE ANGLES* DEGREES 
3.0088E+02 3.0088E+02 

ROTOR SLOPE VECTORS 

J. 2723E-04 2.3270E-14 I.2723E-04 

ROTOR SLOPE VECTOR PHASE ANGLE* DEGREES 

2.2472E+02 2.4623E-25 4.4719E+01 

BEARING SLOPE VECTORS 

0* OOOOE+OO 0. OOOOE+OO 

BEARING SLOPE PHASE ANGLE * DEGREES 
4.5000E+01 4. 5000E+01 

MOUNT SLOPE VECTORS 

I.2723E-04 1.2723E-04 

MOUNT SLOPE PHASE ANGLE* DEGREES 
2*2472E*92 4»47’9E- k 0! 

BEARING FORCE IN X-DIRECTION 
-5. 1425E+01 -5.1425E+0I 

BEARING FORCE IN Y-DIRECTION 
8.5980E+01 3.5980E+01 

BEARING MOMENT IN X-Z PLANE 
0. OOOOE+OO 0. OOOOE+OO 
BEARING MOMENT IN Y-Z PLANE • 

0. OOOOE+OO 0. OOOOE+OO 

MOUNT FORCE IN X-DIRECTION 
I • 7824E+0 I 1 • 7824E+0 1 

MOUNT FORCE IN Y-DIRECTION 
-2.9801E+01 -2.980IE+01 

MOUNT MOMENT IN X-Z PLANE 
0. OOOOE+OO 0. OOOOE+OO 

MOUNT MOMENT IN Y-Z PLANE 
0. OOOOE+OO 0. OOOOE+OO 

BEARING MASS FORCE IN X-DIRECTION 
6.9249E+01 6.9249E+01 

BEARING MASS FORCE IN Y-DIRECTION 
-1.1S78E+02 - I • 1 578E+02 

BEARING INERTIA MOMENT IN X-Z PLANE 
0. OOOOE+OO 0. OOOOE+OO 

BEARING INERTIA MOMENT IN Y-Z PLANE 
0. OOOOE+OO 0. OOOOE+OO 


Figure 46. Computer Starting Rotor Deflection Coordinates 
Run 3 



THE AVERAGE REAL STEP-TIME FOR THIS PRINTOUT = 2.500E-0S 

REAL TIME = 1 . OOOE-03 • SEC SEC 

REVOLUTIONS ARRAYS 

I • 59 1 SE+OO I.59I5E+00 1.59I5E+00 

SPIN SPEED ARRAY* RPMl 

9.5493E*04 9» 5493E+04 9.S493E+04 

ROTOR DISPLACEMENT VECTOR ARRAY. IN 

1 • 01 67E-04 9.3439E-05 1.0167E-04 

ROTOR VECTOR PHASE ANGLE ARRAY. DEGREES 
7 • 3673E+0 1 7.8074E+01 7»3673E+0I 

BEARING DISPLACEMENT VECTORS 
1.0018E-04 U0018E-04 

BEARING DISPLACEMENT PHASE ANGLES. DEGREES 
6.3839E+01 6.3839E+0I 

MOUNT DISPLACEMENT. VECTOR ARRAY. IN 

i • I J « • i O O (.C 4 v J 

MOUNT VECTOR PHASE ANGLE ARRAY* DEGREES 
US385E+02 1 * 5385E+02 

ROTOR WHIRL/SPIN EREO* RATIO ARRAY 

1 • OOCOE+OO l.OOOOE+OO 1 • OOOOE+OO 

ROTOR SLOPE VECTORS 

1 «272 ! E-04 2# 8000E- 1 4 ’ W2721E-04 

ROTOR SLOPE PHASE ANGLES* DEGREES 

7.7677E+01 9*7561E+01 2.5768E+02 

BEARING SLOPE VECTORS 

2.2888E-I5 1.9860E-15 

BEARING SLOPE PHASE ANGLES* DEGREES 
2* 0283E+02 3. 3343E+02 

MOUNT SLOPE VECTORS 

I.2721E-04 N2721E-04 

MOUNT SLOPE VECTOR PHASE ANGLES* DEGREES 
7.7677E+01 2* 5768E+02 


Figure 46 . Concluded 



THE COMPUTED STARTING ROTOR DEFLECTION COORDINATES ARE 


VECTOR ARRAY, IN: 

1.1050E-04 9 • 323 t E- 05 1.1050E-04 

PHASE ANGLE ARRAY, DEGREES: 

2.3I73E+02 2.2480E+02 2.3173E+02 

BEARING DISPLACEMENT VECTOR, IN 
I • 3083E-04 1.3033E-04 

BEARING DISPLACEMENT PHASE ANGLE, DEGREES 
2.4357E+02 2* 4357E+02 

MOUNT VECTORS, IN 

3.2064E-05 3.2064E-05 

MOUNT VECTOR PHASE ANGLES, DEGREES 
1.08S7E+02 1.0857E+02 

R0T0R SLOPE VECTORS 

1.2791E-04 8.6063E-14 I.279IE-04 

ROTOR SLOPE VECTOR PHASE ANGLE, DEGREES 
2.2S48E+02 2.3726E+02 4.S476E+0I 

BEARING SLOPE VECTORS 

O.OOOOEfOO O.OOOOEfOO 
BEARING SLOPE PHASE ANGLE , DEGREES 
4.5000E+01 4. 5000E+0 l 

MOUNT SLOPE VECTORS 

1*2791 E- 04 1.2791E-04 

MOUNT SLOPE PHASE ANGLE, DEGREES 
2.2548E+02 4.S476E+01 

BEARING FORCE IN X-DIRECTI0N 
S.8926E+01 5. 8926E+0 1 

BEARING FORCE IN Y-DIRECTION 
•1*7 538E+02 -1.7538E+02 

BEARING MOMENT IN X-Z PLANE 
0. 0000E+00 O.OOOOEfOO 

BEARING MOMENT IN Y-Z PLANE 
0* OOOOE+OO O.OOOOEFOO • 

MOUNT FORCE IN X-DIRECTION 
-2.0424E+01 -2.0424E+01 

MOUNT FORCE IN Y-DIRECTION 
6* 0789E+0 1 6.0789E+01 

MOUNT MOMENT IN X-Z PLANE 

0* OOOOE+OO O.OOOOEfOO 

MOUNT MOMENT IN Y-Z PLANE 

O.OOOOEfOO O.OOOOEFOO 

BEARING MASS FORCE IN X-DIRECTION 
-7 « 9350E+0 1 -7.93SOEF01 

BEARING MASS FORCE IN Y-DIRECTION 
2.3617EF02 2.361 7E+02 

BEARING INERTIA MOMENT IN X-Z PLANE 
O.OOOOEfOO O.OOOOEFOO 

BEARING INERTIA MOMENT IN Y-Z PLANE 
O.OOOOEfOO O.OOOOEfOO 


i 

Figure 47. Computed Starting Rotor Deflection Coordinates 
Run 4 



THE AVERAGE REAL STEP-TIME FOR THIS PRINTOUT s 2.500E-0S 
REAL TIME * 1.000E-03 SEC SEC 

REVOLUTIONS ARRAY J 

1 • 59 1 5E + 00 l.S915E«-00 1.5915E+00 

SPIN SPEED ARRAY* RPMt 

9» 5493E*04 9.S493E+04 9.5493E+04 

ROTOR DISPLACEMENT VECTOR ARRAY* IN 

I • 1050E-04 9 • 3233E- 05 1.J050E-04 

ROTOR VECTOR PHASE ANGLE ARRAY* DEGREES 
8* 4689E+0 1 7 • 7762E+0 1 8.4689E+0I 

BEARING DISPLACEMENT VECTORS 
I • 3082E-04 1.3082E-04 

BEARING DISPLACEMENT PHASE ANGLES* DEGREES 
9.6528E+0I 9.6528E+01 

MOUNT DISPLACEMENT VECTOR ARRAY* IN 
3.206IE-05 3 • 206 1 E-OS 

MOUNT VECTOR PHASE ANbLE ARRAY* DEGREES 
3.2153E+02 . 3.2153E+02 

ROTOR WHIRL/SPIN FREQ. RATIO ARRAY 

I • OOOOE+OO I » OOOOE+OO I. OOOOE+OO 

ROTOR SLOPE VECTORS 

1 • 2790E-04 I • 1 185E-13 1.2790E-04 

ROTOR SLOPE PHASE . ANGLES* DEGREES 

7.8434E+01 3* 5890E+02 2.5843E+02 

BEARING SLOPE VECTORS 

3.2196E-15 3.3307E- 15 

BEARING SLOPE PHASE ANGLES* DEGREES 
1 • 8000E+02 1-7202E-24 

MOUNT SLOPE VECTORS 

I • 2790E-04 I.2790E-04 

MOUNT SLOPE VECTOR PHASE ANGLES* DEGREES 
7 .8434E+01 2.S843E+02 


Figure 47 . Concluded 



the computed starting rotor deflection coordinates are 


VECTOR ARRAY* IN: 

I.3581E-04 9.2S85E-0S 1.358IE-04 

PHASE ANGLE ARRAY* DEGREES: 

2.2500E+02 2 • 2500E+02 2.2500E+02 

BEARING DISPLACEMENT VECTOR* IN 
2.0785E-04 2.0785E-04 

BEARING DISPLACEMENT PHASE ANGLE* DEGREES 
2.2500E+02 2.2500E+02 

MOUNT VECTORS* IN 

7.204IE-05 7.2041E-05 

MOUNT VECTOR PHASE ANGLES* DEGREES 
4.5000E+01 4.5000E+0I 

ROTOR SLOPE VECTORS 

1 • 3005E-04 2.0684E-14 1.30Q5E-04 

ROTOR SLOPE VECTOR PHASE ANGLE* DEGREES 
2 • 2500E+02 9.0000E+01 4.5000E+01 

BEARING SLOPE VECTORS 

O.OOOOE+OO O.OOOOE+OO 
BEARING SLOPE PHASE ANGLE » DEGREES 
4.S000E+01 4.5000E+01 

MOUNT SLOPE VECTORS 

1 • 3005E-04 1 « 3005E-.04 

MOUNT SLOPE PHASE ANGLE, DEGREES 
2.2500E+02 4.5000E+01 

BEARING FORCE IN X-DIRECTION 
-2.9394E+02 -2.9394E+02 

BEARING FORCE IN Y-DIRECTION 
-2.9394E+02 -2.9394E+02 

BEARING MOMENT IN X-Z PLANE 
O.OOOOE+OO O.OOOOE+OO 

BEARING MOMENT IN Y-Z PLANE 
O.OOOOE+OO O.OOOOE+OO . 

MOUNT FORCE IN X-DIRECTION 
1.0I88E+02 I • 0 1 88E+02 

MOUNT FORCE IN Y-DIRECTION 
1.0188E+02 1 . 0 l 88E+02 

MOUNT MOMENT IN X-Z PLANE 

O.OOOOE+OO O.OOOOE+OO 

MOUNT MOMENT IN Y-Z PLANE 

O.OOOOE+OO O.OOOOE+OO 

BEARING MASS FORCE IN X-DIRECTION 
3. 9582E+02 3.9582E+02 

BEARING MASS FORCE IN Y-DIRECTION 
3.9582E+02 3.9582E+02 

BEARING INERTIA MOMENT IN X-Z PLANE 
O.OOOOE+OO O.OOOOE+OO 

BEARING INERTIA MOMENT IN Y-Z PLANE 
O.OOOOE+OO O.OOOOE+OO 


Figure 48. Computed Starting Rotor Deflection Coordinates 
Run 5 



THE AVERAGE REAL STEP-TIME FOR THIS PRINTOUT * 2.500E-05 

REAL TIME = l • 0O0E-03 SEC SEC 

REVOLUTIONS ARRAY: 

1 • 59 1 5E+00 , 1 • 59 1 5E+00 I.591SE+00 

SPIN SPEED ARRAY* RPM: 

9.5493E+04 9.5493E+04 9.5493E+04 

ROTOR DISPLACEMENT VECTOR ARRAY* IN 

1 » 3580E-04 9.2587E-0S 1.3580E-04 

ROTOR VECTOR PHASE ANGLE ARRAY* DEGREES 
7 • 7958E+0 1 7.79S9E+01 7.7958E+01 

BEARING DISPLACEMENT VECTORS 
2» 0783E- 04 2.0783E-04 

BEARING DISPLACEMENT PHASE ANGLES* DEGREES 
7 • 7957E+0 1 7.7957E+01 

MOUNT DISPLACEMENT VECTOR ARRAY* IN 
7 « 2C27E-C5 7.SC27S-05 

MOUNT VtCTUR" PHASE ANGLE ARRAY* DEGREES 
2.5796E+02 2.5796E+02 

ROTOR WHIRL/SPIN FREO. RATIO ARRAY 

1 • OOOOE+OO 1 • OOOOE+OO l.OOOOE+OO 

ROTOR SLOPE VECTORS 

l*3003E-04 2.7882E-I4 1.3003E-04 

ROTOR SLOPE PHASE ANGLES* DEGREES 

7.7951E+01 1 • 9796E+02 2. S796E+02 

BEARING SLOPE VECTORS 

S. J837E-15 4.0943E-15 .■*' ; . 

BEARING SLOPE PHASE ANGLES* DEGREES 
2* 2326E+02 3.1940E+02 

MOUNT SLOPE VECTORS 

1.3003E-04 1 • 3003E-04 

MOUNT SLOPE VECTOR PHASE ANGLES* DEGREES 
7.7957E+01 2.5796E>02 


Figure 48. Concluded 



THE COMPUTED STARTING ROTOR DEFLECTION COORDINATES ARE: 
VECTOR ARRAY* IN: 

1.1047E-04 9.041SE-05 1.1047E-04 

PHASE ANCLE ARRAY* DEGREES: 

2.2500E+02 2.2500E+02 2.2500E+02 

BEARING DISPLACEMENT VECTOR* IN 
7 • 3652E- 05 7.3652E-05 

BEARING DISPLACEMENT PHASE ANGLE* OEGREES 
2.2500E+02 2.2500E+G2 

MOUNT VECTORS* IN 

3.6S26E-0S 3.6R26E-05 

MOUNT VECTOR PHASE ANGLES* DEGREES 
2.2500E+02 2.2S00E+02 

ROTOR SLOPE VECTORS 

I • 8239E-04 1.1243E-14 1.S239E-04 

ROTOR SLOPE VECTOR PHASE ANGLE* DEGREES 
2.2500E+02 2.6069E+02 4.5000E+01 

BEARING SLOPE VECTORS 

2.2062E-04 2.2062E-04 

BEARING SLOPE PHASE ANGLE * DEGREES 
2.2500E+02 4. 5000E+0 I 

MOUNT SLOPE VECTORS 

3.8234E-05 3.8234E-05 

MOUNT SLOPE PHASE ANGLE* DEGREES 
4 » 50CCE+0 1 2.C5CCC+02 

BEARING FORCE IN X-DIRECTION 
-5.2080E+01 -5.2080E+01 

BEARING FORCE IN Y-DIRECTION 
-5.2080E+01 -5.2080E+01 

BEARING MOMENT IN X-Z PLANE 
- I « 5600E+02 1 • 5600E+02 

BEARING MOMENT IN Y-Z PLANE 
-1.5600E+02 1.5600E+02 

MOUNT FORCE IN X-DIRECTION 
-5.2080E+0I -5.2080E+01 

MOUNT FORCE IN Y-DIRECTION 
-5.2080E+01 -5.2080E+01 

MOUNT MOMENT IN X-Z PLANE 
5.4072E+01 -5.4072E+01 

MOUNT MOMENT IN Y-Z PLANE 
5.4072E+0I -5* 4072E+0 I 

BEARING MASS FORCE IN X-DIRECTION 
O.OOOOE+OO O.OOOOE+OO 

BEARING MASS FORCE IN Y-DIRECTION 
O.OOOOE+OO O.OOOOE+OO 

BEARING INERTIA MOMENT IN X-Z PLANE 
2. I008E+02 -2.I008E+02 

BEARING INERTIA MOMENT IN Y-Z PLANE 
2. 1 008E+02 -2. I 008E+02 


Figure 49. Computed Starting Rotor Deflection Coordinates 
Run 6 



THE AVERAGE REAL STEP-TIME F0R THIS PRINTOUT * 2.SOOE-OS 
REAL TIME = I.O00E-03 SEC SEC 

REVOLUTIONS ARRAY: 

1*591 SE+OO 1.S915E+00 1.5915E+00 

SPIN SPEED ARRAY* RPM: 

9.5493E+04 9.5493E*04 9.5493E+04 

ROTOR DISPLACEMENT VECTOR ARRAY* IN 

1*1 047E-04 9.0420E-05 1.1047E-04 

ROTOR VECTOR PHASE ANGLE ARRAY* DEGREES 
7 • 7958E+0 1 7.7959E+01 7.7958E+01 

BEARING DISPLACEMENT VECTORS 
7*3651E-05 7.3651E-05 

BEARING DISPLACEMENT PHASE ANGLES* DEGREES 
7 • 7958E+G 1 7.7958E+01 

MOUNT DISPLACEMENT. VECTOR .ARRAY* IN 
3 j 6826E-C5 C.632yE-05 

MOUNT VECTOR PHASE ANGLE ARRAY* DEGREES 
7 • 7958E+0 1 7.7958E+01 

ROTOR WHIRL/SPIN FRES* RATIO ARRAY 

1 •OOOOE+OO 1 • OOOOE+OO l.OOOOE+OO 

ROTOR SLOPE VECTORS 

1 • 8236E-04 7.9395E-15 1.8236E-04 

ROTOR SLOPE PHASE ANGLES* DEGREES 

7.7952E+OI 4.5S79E+01 2. 5795E+02 

BEARING SLOPE VECTORS 

2.2058E-04 2.2058E-04 

BEARING SLOPE PHASE ANGLES* DEGREES 
7.7951E+01 2* S79 5E+02 

MOUNT SLOPE VECTORS . r" ’ 

3.8224E-05 3.8224E-05 

MOUNT SLOPE VECTOR PHASE ANGLES* DEGREES 
2.5795E+02 7.7946E-MH 


Figure 49. Concluded 



Tie computed starting rotor deflection coordinates are? 

VECTOR ARRAY# IN: 

I.0867E-04 9.3S46E-05 t.0867E-04 

PHASE ANGLE ARRAY# DEGREES t 

2.2461E+02 2.258SE+02 2.2461E+02 

BEARING DISPLACEMENT VECTOR# IN 
7.2453E-05 7 • 2453E-05 

BEARING DISPLACEMENT PHASE ANGLE# DEGREES 
2. 2461 E + 02 2.2461E+02 

MOUNT VECTORS# IN 

3.6226E-05 3.6226E-05 

MOUNT VECTOR PHASE ANGLES# DEGREES 
2.2461E+02 2. 2461 E+02 

ROTOR SLOPE VECTORS 

I • I 960E-04 1.2668E-J4 1.1960E-04 

ROTOR SLOPE VECTOR PHASE ANGLE# DEGREES 
2.1202E+02 3.2791E+02 3.2025E+01 

BEARING SLOPE VECTORS 

I • 1 785E-04 1 • I 785E-04 

BEARING SLOPE PHASE ANGLE # DEGREES 
2.0219E+02 2.2193E+0I 

MOUNT SLOPE VECTORS 

2.C424r-OS 2-0424E-05 

“SUNT SLOPE PHASE ANGLE# DEGREES 
2.9219E+02 i • 121 9E+02 

BEARING FORCE IN' X-DIRECTION 
-5.1577E+01 -5. 1577E+0J 

BEARING FORCE IN Y-DIRECTION 
-5.0884E+01 -S.0884E+01 

BEARING MOMENT IN X-Z PLANE 
-4. 4515E+01 4.4515E+01 

BEARING MOMENT IN Y-Z PLANE 
1 • 09 12E+02 -1.0912E+02 

MOUNT FORCE IN X-DIRECTION 
-5.IS77E+01 -5.I577E+01 

MOUNT FORCE IN Y-DIRECTION 
-5.0884E+01 -5.0884E+01 

MOUNT MOMENT IN X-Z PLANE 
1.5429E+01 - 1 . 5429E+0 I 

MOUNT MOMENT IN Y-Z PLANE 

-3.7822E+01 3.7822E+0! 

BEARING MASS FORCE IN X-DIRECTION 
O.OOOOE+OO O.OOOOE+OO 

BEARING MASS FORCE IN Y-DIRECTION 
O.OOOOE+OO O.OOOOE+OO 

BEARING INERTIA MOMENT IN X-Z PLANE 
S.994SE+01 -S.994SE+0I 

BEARING INERTIA MOMENT IN Y-Z PLANE 
- t » 4694E+02 1 c 4694E+02 


Figure 50. Computed Starting Rotor Deflection Coordinates 
Run 7 
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THE AVERAGE REAL STEP-TIME FOR THIS PRINTOUT = 2.500E-05 

SEC 

REAL TIME » I.000E-03 SEC 

REVOLUTIONS ARRAY* 

1.S915E+00 l • 59 1 5E*00 1.5915E+00 

SPIN SPEED ARRAY* RPMl 

9 • 5493E+04 9.S493E+04 9.S493E+04 

ROTOR DISPLACEMENT VECTOR ARRAY* IN 

1.0867E-04 9.3848E-05 1.0867E-04 

ROTOR VECTOR PHASE ANGLE ARRAY* DEGREES 
7.7571E+01 7.8809E+0I 7.7S71E+01 

BEARING DISPLACEMENT VECTORS 
7 • 2452E-05 7.2452E-05 

BEARING DISPLACEMENT PHASE ANGLES* OEGREES 
7.7571E+01 7 • 7 5 7 1 E* 0 1 

MOUNT DISPLACEMENT VECTOR ARRAY* IN 
3.6226E-05 3.6226E-05 

MOUNT VECTOR PHASE ANGLE ARRAY* DEGREES 
7.7571E+01 7.7571E+01 

ROTOR WHIRL/SPIN FREQ. RATIO ARRAY 

I » OOOOE+OO 9.9999E-0I t.OOOOE+OO 

ROTOR SLOPE VECTORS 

I • I 9 59 E- 04 I.7485E-14 I.1959E-04 

ROTOR SLOPE PHASE ANGLES* DEGREES 

6.4989E+01 9.26S7E+01 2.4499E+02 

BEARING SLOPE VECTORS 

1 • 1 783E-04 1 • 1 783E-04 

BEARING SLOPE PHASE ANGLES* DEGREES 
5*5161 E+ 01 2.3516E+02 

MOUNT SLOPE VECTORS 

2.0412E-05 2.0412E-05 f >i.- 

MOUNT SLOPE VECTOR PHASE ANGLES* DEGREES 
1 *451 5E+02 3.2515E»02 


Figure 50. Concluded 



THE COMPUTED STARTING R0T0R DEFLECT ION COORDINATES ARE 
VECTOR ARRAY, INt 

I.0699E-04 9.7030E-05 1.0699E-04 

PHASE ANGLE ARRAY, DEGREES: 

2.2526E +02 2.2446E+02 2.2526E+02 

BEARING DISPLACEMENT VECTOR, IN 
7.1331E-05 7*1331 E-05 

BEARING DISPLACEMENT PHASE ANGLE, DEGREES 
2* 2526E+02 2.2526E+02 

MOUNT VECTORS, IN 

3.5665E-05 3.5665E-05 

MOUNT VECTOR PHASE ANGLES, DEGREES 
2* 2526E+02 2.2526E+02 

ROTOR SLOPE VECTORS 

5* 783 1 E- 05 2. 1740E- 14 5.7831E-05 

ROTOR SLOPE VECTOR PHASE ANGLE, DEGREES 
2.4266E+02 2.5045E+02 6.2656E+01 

BEARING SLOPE VECTORS 

3*011 8E-0S 3*011 8E- 05 

BEARING SLOPE PHASE ANGLE , DEGREES 
3*07 1 5E+02 1.2715E+02 

MOUNT SLOPE VECTORS 

5* 2456E-05 5.2456E-05 

MOUNT SLOPE PHASE ANGLE, DEGREES 
2. 1 » 44£+o? 3* 1 443E+0 1 

BEARING FORCE IN X-DIRECTION 
-5.0212E+01 -5*021 2E+0 1 

BEARING FORCE IN Y-DIRECTION 
-5.0664E+01 -5.0664E+01 

BEARING MOMENT IN X-Z PLANE 
2.S824E+02 -2.5824E+02. 

BEARING MOMENT IN Y-Z PLANE 
I • 5790E+02 “ 1 • 5790E+02 

MOUNT FORCE IN X-DIRECTION 
-5.0212E+01 -5* 02 12E+0I 

MOUNT FORCE IN Y-DIRECTION 
- 5* 0664E+0 1 - 5* 0664E + 0 1 

MOUNT MOMENT IN X-Z PLANE 
-8.9507E+0I 8.9507E+01 

MOUNT MOMENT IN Y-Z PLANE 
-5.4728E+0I 5.4728E+01 

BEARING MASS FORCE IN X-DIRECTION 
O.OOOOE+OO O.OOOOE+OO 

BEARING MASS FORCE IN Y-DIRECTION 
O.OOOOE+OO O.OOOOE+OO 

BEARING INERTIA MOMENT IN X-Z PLANE 
-3.4775E+02 3* 4775E+02 

BEARING INERTIA MOMENT IN Y-Z PLANE 
-2* I262E+02 2* 1262E+02 WHAT? 


Figure 51. Computed Starting Rotor Deflection Coordinates 
Run 8 



THE AVERAGE REAL STEP-TIME FOR THIS PRINTOUT = 2.500E- 

REAL TIME = 1-000E-03 SEC 05 SEC 

REVOLUTIONS ARRAYl 

1.5915E+00 1 • 59 1 5E+00 I.5915E+00 

SPIN SPEED ARRAY* RPMl 

9.5493E+04 9*5493E*04 9.5493E+04 

ROTOR DISPLACEMENT VECTOR ARRAY* IN 

I.0699E-04 9 • 703 1 E- 05 U0699E-04 

ROTOR VECTOR PHASE ANGLE ARRAY* DEGREES 
7.8215E+01 7.7421E+01 7.8215E+01 

BEARING DISPLACEMENT VECTORS 
7.1330E-0S 7 » 1330E-05 

BEARING DISPLACEMENT PHASE ANGLES* DEGREES 
7.8215E+01 7.8215E+0I 

..MOUNT DISPLACEMENT VECTOR ARRAY* IN 
3.5665E-05 . 3.5665E-05 

MOUNT VECTOR PHASE ANGLE ARRAY* DECREES 
7*821 5E+01 7*821 5E+0 l 

ROTOR WHIRL/SPIN FREO* RATIO ARRAY 

1 *OOOOE+O0 1 • OOOOE+OO l.OOOOE+OO 

ROTOR SLOPE VECTORS 

5.7828E-05 2.7954E-14 5.7828E-0S 

ROTOR SLOPE PHASE ANGLES* DEGREES 

9*561 6E+01 3.5497E+02 2.7S62E+02 

BEARING SLOPE VECTORS 

3.0120E-05 3.0120E-05 

BEARING SLOPE PHASE ANGLES* DEGREES 
1.601 1E+02 3*401 1 E+02 

MOUNT SLOPE VECTORS 

5.2454E-05 5.2454E-05 

MOUNT SLOPE VECTOR PHASE ANGLES* DEGREES 
6* 4400E+0I 2. 4440E+02 


Figure 51; Concluded 



the computed starting rotor deelection coordinates are: 

VECTOR ARRAY* IN: 

1.0636E-04 9.8228E-0S 1.0636E-04 

PHASE ANGLE ARRAY* DEGREES: 

2.2500E+02 2. 2500E+ 02 2.2SOOE+02 

BEARING DISPLACEMENT VECTOR* IN 
7 • 0903E-05 7.0908E-05 

BEARING DISPLACEMENT PHASE ANGLE* DEGREES 
2.2S00E+02 2. 2500E+02 

MOUNT VECTORS* IN 

3.54S4E-05 3 • 5454E- 05 

MOUNT VECTOR PHASE ANGLES* DEGREES 
2. 2500E+02 2.2500E+02 

ROTOR SLOPE VECTORS 

3* 1 959E-05 1.1607E-14 3.1959E-05 

ROTOR SLOPE VECTOR PHASE ANGLE* DEGREES 
2.2500E+02 3 • 0220E+02 4.5000E+01 

BEARING SLOPE VECTORS 

3.S262E-05 3.5262E-05 

BEARING SLOPE PHASE ANGLE * DEGREES 
4.5000E+0J 2.2500E+02 

MOUNT SLOPE VECTORS 

6.7222'E-OS 6.7222E-05 

MOUNT SLOPE PHASE ANGLE* DEGREES 
2.2500E+02 4. 5000E+01 

BEARING FORCE IN X-DIRECTION 
-S.0139E+01 -S.0139E+0J 

BEARING FORCE IN Y-DIRECTI0N 
-S.0I39E+0I -5.0139E+0I 

BEARING MOMENT IN X-Z PLANE 
2.7428E+02 -2.7428E+02 

BEARING MOMENT IN Y-Z PLANE 
2.7428E+02 -2.7428E+02 

MOUNT FORCE IN X-DIRECTION 
-S.0139E+01 -5.0139E+01 

MOUNT FORCE IN Y-DIRECTION 
-S.0139E+0I -5.0139E+0I 

MOUNT MOMENT IN X-Z PLANE 
-9.5066E+0I 9.5066E+01 

MOUNT MOMENT IN Y-Z PLANE 
-9. 5066E+01 9 • S066E+0 1 

BEARING MASS FORCE IN X-DIRECTION 
O.OOOOE+OO O.OOOOE+OO 
BEARING MASS FORCE IN Y-DIRECTION 
0* OOOOE+OO O.OOOOE+OO 

BEARING INERTIA MOMENT IN X-Z PLANE 
-3* 6934E+02 3. 6934E+02 

BEARING INERTIA MOMENT IN Y-Z PLANE 
-3.6934E+02 3.6934E+02 


Figure 52. Computed Starting Rotor Deflection Coordinates 
Run 9 



THE AVERAGE REAL STEP-TIME FOR THIS PRINTOUT a 2.500E-05 SE 
REAL TIME * l • 000E-03 SEC 
REVOLUTIONS ARRAY I 

1.591SE+00 1.5915E+00 1.5915E+00 

SPIN SPEED ARRAY* RPM» 

9» 5493E+04 9.5493E+04 9.5493E-MJ4 

ROTOR DISPLACEMENT VECTOR ARRAY# IN 

1 • 0636E-04 9.8228E-05 I.0636E-04 

ROTOR VECTOR PHASE ANGLE ARRAY# DEGREES 
7.7958E+01 7.7958E+01 7.7958E+0I 

BEARING DISPLACEMENT VECTORS 
7.0907E-05 7.0907E-05 

BEARING DISPLACEMENT PHASE ANGLES* DEGREES 
7.7958E+0I 7.7958E+0I 

MOUNT DISPLACEMENT VECTOR ARRAY# IN 
3.S454E-05 3.5454E-05 

MOUNT VECTOR PHASE ANGLE ARRAY* DEGREES 
7.7958E*01 7 • 7958E+0 1 

ROTOR WHIRL/SPIN FREO. RATIO ARRAY 

t*OOOOE+QO 1 •OOOQE+OO l.OOOOE+OO 

ROTOR SLOPE VECTORS 

3.1955E-05 I • 5493E- 1 4 3.I95SE-05 

ROTOR SLOPE PHASE ANGLES# DEGREES 

7.796SE+01 4» 3 1 44E+0 1 2.5796E+02 

BEARING SLOPE VECTORS 

3.5267E-05 3.5267E-05 

BEARING SLOPE PHASE ANGLES* DEGREES ; 

2.5795E+02 7.7948E+0I 

MOUNT SLOPE VECTORS 

6* 7222E-05 6.7222E-05 

MOUNT SLOPE VECTOR PHASE ANGLES* DEGREES 
7.7956E+0! 2.S796E+02 


Figure 52. Concluded 
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THE COMPUTED STARTING ROTOR DEFLECTION COORDINATES ARE: 


VECTOR ARRAY* INI 

8.8108E-05 9.9668E-05 8.8108E-05 

PHASE ANGLE ARRAY* DEGREES: 

2.6033E+02 2.2282E+02 2. 6033E+02 

BEARING DISPLACEMENT VECTOR* IN 
5*161 9E- 05 5*161 9E- 05 

BEARING DISPLACEMENT PHASE ANGLE* DEGREES 
3*057 1E+02 3.0571E+-02 

MOUNT VECTORS* IN 

6.3544E-05 6.3544E-05 

MOUNT VECTOR PHASE ANGLES* DEGREES 
2.2501E+02 2.2501E+02 

ROTOR SLOPE VECTORS 

5* 6783E- 05 1.3875E-14 5.6783E-0S 

ROTOR SLOPE VECTOR PHASE ANGLE* DEGREES 
2* 9533E+02 2.8458E+02 1.1533E+02 

BEARING SLOPE VECTORS 

3.3267E-05 3.3267E-05 

BEARING SLOPE PHASE ANGLE * DEGREES 
3.4071E+02 !.6071P>02 

MOUNT SLOPE VECTORS 

4.0953E-05 4.0953E-05 

MOUNT SLOPE PHASE ANGLE* DEGREES 
2* 6001 E+02 8.Q012E+01 

BEARING FORCE IN X-DI RECTION 
7.0864E+02 7.0864E+02 

BEARING FORCE IN Y-DIRECTION 
-1.8988E+02 - 1 .8988E+-02 

BEARING MOMENT IN X-Z PLANE 
4.4431E+02 -4.4431E+02 

BEARING MOMENT IN Y-Z PLANE 
1 • 6 1 69E+02 -1.6169E+02 

MOUNT FORCE IN X-DI RECTION 
3* 5959E+02 3.5959E+02 

MOUNT FORCE IN Y-DIRECTION 
-5*391 IE+02 -5*3911 E+02 

MOUNT MOMENT IN X-Z PLANE 
3.8912E+02 -3*89 12E+02 

MOUNT MOMENT IN Y-Z PLANE 
- I.S170E+02 1*51 70E+02 

BEARING MASS FORCE IN X-DIRECTION 
-3» 4905E+02 -3.4905E+02 

BEARING MASS FORCE IN Y-DIRECTION 
-3.4922E*02 -3 * 4922E+02 

BEARING INERTIA MOMENT IN X-Z PLANE 
-5* 5194E+01 5. 5194E*'01 

BEARING INERTIA MOMENT IN Y-Z PLANE 
-.3* 1 339E+02 3. 1339E+02 


Figure 53. Computed Starting Rotor Deflection Coordinates 
Run 10 



THE AVERAGE REAL STEP-TIME FOR THIS PRINTOUT *2. 
REAL TIME = 4.000E-03 SEC 

REVOLUTIONS ARRAY J 

6.3662E+00 6.3662E+00 6.3662E+00 

SPIN SPEED ARRAY* RPM* 

9.5493E+04 9.5493E+04 9.5493E+04' 


ROTOR DISPLACEMENT VECTOR ARRAY* IN 

8* 8 108E-05 9.9668E-05 8.8108E-05 

ROTOR VECTOR PHASE ANGLE ARRAY* DEGREES 
3.2I67E+01 3.546SE+02 3.2167E+01 

BEARING DISPLACEMENT VECTORS 
S* ! 623E-05 5.I623E-05 

BEARING DISPLACEMENT PHASE ANGLES* DEGREES 
7.7545E+01 7.7545E+OI 

MOUNT DISPLACEMENT VECTOR ARRAY* IN 
6.3547E-0S 6.3S47E.-05 

i *r/» n «-»* i « o i- a r~ a no a v r\ r* r* o r» r* f" 

Wiv I iuh *-.iy uu L. nhi.nlJ vuU(\uu J 

3* 3684E+02 3.S6B4E+02 


ROTOR WHIRL/SPIN FREQ. RATIO ARRAY 

9.999SE-01 ! • OOOOE+OO 9.9995E-01 . 

BEARING .MASS WHIRL/ROTOR SPIN FREQ. RATIO ARRAY* 
1.0001E+00 1.0001E+00 

ROTOR SLOPE VECTORS 

5.678SE-05 1.2678E-14 5.678SE-0S 

ROTOR SLOPE PHASE ANGLES* DEGREES 

6* 7 I 62E+0 1 3. 4538E+02 2.4716E+02 

BEARING SLOPE VECTORS 

3.3270E-05 3.3270E-05 ^ 

BEARING SLOPE PHASE ANGLES* DEGREES 
! . 1253E+02 2.9254E+02 

MOUNT SLOPE VECTORS 

4.0954E-05 4.0954E-0S 

MOUNT SLOPE VECTOR PHASE ANGLES* DEGREES 
3. I837E+01 2. I 184E+02 


Figure 53. Concluded 
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STATION NO : 1 


2 


3 



YOUNG'S MODULUS . - 20.68xl0 6 NENT0NS/C2^(3OxlO 6 pai) 

SHEAR MODULUS 8.27x10° NEW0NS/CM 2 (l2.xl0 6 psi) 

POISSON'S RATIO =0.3 

MATERIAL DENSITY - 0 

ROTOR MASS = 22.68 Kg (5Q LB)/STATION 

ROTOR TRANSVERSE MASS INERTIA = 146.3 Kg-qf (50 LB-IN 2 )/STATION 
ROTOR POLAR MASS INERTIA = 87.8 Kg-CM 2 (30 LB-IN 2 )/STATI0N 

ROTOR ECCENTRICITY =0.000254 CM/STATION 

ECCENTRICITY PHASE ANGLE . = 45° /STATION. 

ROTOR INERTIA INERTIA MISALIGNMENT =0.0057°, 0°, -0.0057° 
MISALIGNMENT PHASE ANGLE = 45°/STATI0N' 

ROTOR SPIN SPEED = 95493. RPM 


Figure 53. Schematic Representation, of the Rotor System and 
Properties Used in the Computer Runs 




IV. FINAL VERIFICATION OF IBM 360/370 COMPUTER PROGRAM 


After conversion of the GE computer program to an IBM 360/370 version, the 
IBM program was checkout in five different runs: 

1. Steady-state run in English units 

2. Transient spin speed run in English units 

3. Steady-state run in international units 

4. Transient spin speed run in international units (Appendix F, 

Table XXIII) 

5. A 15-rotor station and 6 bearing station run using the current maxi- 
mum computer program capacity (Appendix F, Table XXIV) 

Runs from 1 through 4 were also made on the GE computer version using identical 
inputs. Basically identical results, except round off errors, were observed 
for both computer versions. 

The same physical rotor design and operation data were used in the computer 
runs involving English and international units. Equivalent computation re- 
sults were obtained for the two different systems of units. The IBM computer 
results for the final computer checkout using Runs 4 and 5 are attached as 
Appendix F (Tables XXIII and XXIV, respectively). In all the IBM 360/370 runs, 
the Runge-Kutta integration technique was used instead of the Adams-Moulton 
technique. For reasonably smooth rotor acceleration function, the former tech- 
nique is faster than the latter. In the case where the rotor acceleration 
function is unpredictable, the Adams-Moulton predictor-corrector technique 
will be applied. 
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CONCLUSIONS AND RECOMMENDATIONS 


During the contract period, extensive efforts to update and broaden the capa- 
bility of the computer program have been made. An optimum integration tech- 
nique and solution method have been established and included in the program. 

The computational speed over that from the previous contractual study has been 
substantially improved by a factor of 2970. The CPU to real-time ratio for a 
15-station and 6-bearing program is 5050. The input/output time versus real- 
time ratio will vary according to the frequency of input/output for a computer 
run. For this run, the input/output (channel) to real-time ratio is 5123. 
Several useful rotor dynamics parameters such as hysteresis, general in-phase 
and out-of-phase bearing force and moment characteristics, transverse effects 
of torsional and axial loading, and bearing mass and inertia effects have been 
included in the program. To facilitate computation processing, flexibility in 
optionally selecting some of the rotor parameters was provided. Program input/ 
output may be performed in international or English units system as a user's 
option. Although the program does not include all possible rotor dynamics 
parameters at present,' it should cover most rotor designs where the applica- 
tion requires stringent design and stable operation. 

Future efforts in experimental verification and updating of the computer pro- 
gram will be of great interest in keeping pace with the demands in rotor dy- 
namics technology. Experimental effort to verify the various applications of 
the transient rotor dynamics analysis computer program is recommended as a 
follow-on to the analytical effort completed in this contract. Periodic up- 
dating of the transient rotor dynamics program should be under taken to keep 
pace with the future technologocal demands and advances in the field of 
rotor dynamics. Accordingly, specific recommendations are delineated below: 

1. Experimental verification of the validity and accuracy of the various 
simulation options in the computer model. This could include rotor 
hysteresis effects simulation. 

2. Rotor casing mass and mass moment of inertia. The. rotor casing can 
be rigid or flexible, rotating or nonrotating. There can be several 
rotor casings with or without coupling between them. A rotating 
casing will in effect be a multiple concentric rotating rotor system. 

3. Rotor casing support. The rotor casing can be supported on a found- 
ation attached to an inertial frame of reference or the foundation 
may experience an angular velocity or acceleration such as space 
vehicle borne rotating machinery. The rotor casing may be floating 
or attached to a floating mass. 

4. Axial elasticity of rotor and resulting axial rotor dynamics. 

5. Compute and writeout rotor dynamic stresses at specified locations. 

6. Design a simplified computer model to include a minimum number of 
basic parameters to be used for quick preliminary exploration 
analysis in a new design or study. 
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APPENDIX A 

COMPUTER PROGRAM USER'S INSTRUCTION 
DESCRIPTION OF THE COMPUTER PROGRAM 

The transient analysis computer program, written in Fortran IV, consists of 
one main program, 11 subroutines, and one Fortran function. Their names and 
the basic flow path are depicted in Fig. 56. ' The functions of the main program, 
the subroutines, and the Fortran function are delineated as follows. 

Main Program 

The main program is used as the basic calling program to coordinate the opera- 
tion of the subroutines and the function so as to perform the various prescribed 
calculations. The other important function performed within the main program 
is to provide English unit output of the computed results in the printout and 
graphical output as required. In addition, it also prints out the following 
nondimens ional output which are common to both English unit and International 
unit outputs: 

1. Average real time- step for this printout 
2., Real time 

3. Rotor spin revolution array. 

Subroutine HYSREA 

The subroutine HYSREA reads all of the input data. It will also -read punched 
cards when C0NTIN = 1 . The namelist read procedure is basically used in 
HYSREA except for (1) a title card which precedes the namelist read statement 
and (2) the punch cards following the namelists. There are two namelists 
defined, MUST and OPTION, with MUST namelist proceeding the OPTION namelist. 
Every Fortran variable except ID is required input in the MUST namelist. ID 
is included only for the purpose of punching card sequence numbers and is not 
used in the computations. The variables included in the OPTION namelist may 
be read at the user's option. For this reason, a set of default values for 
the OPTION namelist, variables are stored in HYSREA which can be overwritten 
by the desired OPTION input data. 

Subroutine HYSWRI 

HYSWRI is a subroutine that prints out the input data including the built-in 
data in HYSREA if it is not overwritten by actual input. The data writeout in 
HYSWRI will include the values of general, nondimensional data and the remain- 
der of the input data when English units are specified. 
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PUNCHED CARDS INPUT 
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Figure 56, Computer Program Components and Flow Path 




Subroutine HYSWME 


HYSWME is another subroutine that prints out the data following the general, 
nondimens ional data which are printed out in HYSWRI. HYSWME is used when input 
and output are required in International units. 

Subroutine HYSMET „ 

HYSMET will write out computed results and produce graphical output when input 
and output are required in International units. 

Subroutine HYSINF 


HYSINF is used to generate rotor displacement and slope deflection influence 
coefficients due to applied unit forces and moments. The influence coefficients 
are referred to a straight line joining the centers of the first and last 
bearings. 

Subroutine HYSSTA 

HYSSTA will compute startup rotor dynamic deflections based on a steady-state 
and axisymmetric rotor-bearing system design. Influence coefficients provided 
by HYSINF will be used in generating the startup rotor configuration through 
the use of the simultaneous equation solution Fortran function. Rotor mass 
and inertial properties are also computed in HYSSTA. 

Subroutines HYSRKA, RUNKUT and ADAMLT 

These are the Adams-Moulton and Runge-Kutta integration subroutines. HYSRKA is 
called by MAIN program and in turn RUNKUT and ADAMLT are called by HYSRKA. 

FUND is called by these subroutines to return with derivatives corresponding to 
the calling variables. There are three optional integration techniques to be 
selected for use according to the value assigned to the Fortran variable, IND, 
the options are: 

IND - 0 Uses Adams-Moulton predictor-corrector variable step integration 
technique. TOLI is used in determining the time steps applied 
for a desired accuracy. 

IND = 1 Uses 4th order Runge-Kutta fixed step integration technique 
IND = 2 Uses Adams-Moulton fixed step integration technique 
Subroutine FUND 

FUND is used to generate time derivatives of the incoming variables from the 
calling program. From the incoming variables, the rotor and bearing reactional 
forces and moments are first determined and the time derivatives are computed 
from the mass properties and the reactional loads. 
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Function ISIMEQ 


ISIMEQ is a simultaneous equation solution function called by HYSSTA in com- 
puting the startup rotor deflection configuration. 

Subroutine PLOTBR 

PLOTBR is called by MAIN or HYSMET for graphical plotting of bearing force and 
displacement versus rotor spin speed function. 

INPUT PROCEDURE 

The computer program is written to simulate a continuous rotor mass distribution 
by a discrete-mass rotor model with an appropriate mass-less elastic shaft. 

In general, the minimum number of discrete masses used should be such that the 
desired rotor dynamic mode shape can be sustained. For instance, if a shaft 
operates in a bending critical speed range, the mode shape could be 
to sustain this mode shape a minimum of 5 masses is required. In general, to 
obtain good accuracy several times the minimum number of mass requirements are 
used. For rotor motion predominately influenced by mass eccentricity and damp- 
ing and stiffness function, the mode shape in a critical speed range may be 
substantially modified from that of a pure critical speed mode shape. Judge- 
ment must hence be exercised in selecting the number of discrete masses to 
adequately represent a rotor configuration. 

The rotor to be studied is first divided into consecutively numbered stations. 
The total number of stations may vary from 3 to 15, inclusively. Rotor sec- 
tions between adjacent rotor stations are labeled with the same numbers as that 
of the left adjacent stations. The rotor property input data are appropriately 
subscripted according to the rotor station or section numbers. For nonlinear 
stiffness bearing data, two-dimensional subscripts are used. The first sub- 
script defines their bearing station location and the second defines the non- 
linear bearing stiffness sections. 

A complete input data writeout is provided for each computer run. A detailed 
description of the input, output, and usage of the program appears in the 
following sections. 


Input 

A namelist input procedure is used as the basic input format due to its flexi- 
bility in selecting input parameters and the liberal use of built-in input data 
when appropriate. For preliminary analysis, by making use of built-in data, 
the input data volume, particularly for a large number of rotor stations, can 
be drastically reduced. 
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The complete input data must include the following sections in the sequence 
listed. 

1. Title. One 80-column card, with the first 72 columns available for a 
descriptive title and columns 73 through 80 reserved for card identi- 
fication, must be provided. 

2. Name list /MUST/. This namelist data section consists of variable names 
which, except for the variable name ID, must be read in. The input 
variables are defined as follows: 

a. Integration step (real) time (DT) , seconds 

b. Maximum run (real) time (TMAX) , seconds 

c. Minimum printing (real) time interval (DP), seconds 

d. Total number of rotor stations (NS) 

e. Total number of bearing stations (NB) 

f . Startup whirl and spin speed (FDOTI) , rpm 

g. Bearing location rotor stations (IB(K)) 

h. Rotor section outside diameters (DD(I)), inches 

i. Rotor section length (QL(I)), inches « 

j. MET = Input/output English or international unit control variable 

k. ID, not a part of rotor dynamics analysis input data. ID is includ 
included for the purpose of allowing data cards to be numbered in 
order to maintain their proper sequence. 

The namelist/MUST/ pertains to rotor geometry and bearing arrangement 
data. It is not possible to provide built-in values to approximate 
the input rotor and bearing information. The physical input of these 
data is necessary. 

3. Namelist/0PTION/. This input data included in this section may be 
applied at the user's option, although a blank card, shown at top of 
next page, is necessary even if no optional input is desired. The 
user needs only to include the data different from built-in values. 

The complete list of namelist/0PTION/ is shown on the following page 
and the data built-in values are described ip subroutine HYSREA in 
Appendix C. 
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4. Punched Card Read In (Conditional Input Data). For a continued analy- 
sis from a previous study punched cards generated from the previous 
analysis must be provided. Concurrently the data variable CONTIN = 1, 
in namelist/OPTION/ must be entered. In addition, the data entered in 
namelist/MUST/ and namelist/OPTION/ must be identical to those from 
the previous study. 

If CONTIN = 0, or no entry of CONTIN is made in namelist/OPTION, 
punched cards must not be included. 

Although the startup rotor dynamic configuration is not used in a 
continued analysis, the startup configuration will still be generated. 
The startup configuration may be used as a verification for the input 
data which must be the same as those for the original run. 

Data names included in namelist/MUST/ and namelist/OPTION/ are as follows. 

NAMELIST/MUST/DT,TMAX,DP,NS,NB,FD0T1, IB,DD > QL,MET, ID 

NAMELIST/OPTION/ IND.TOLI ,T, CONTIN, ITORQ, IPP, IMT, RIG, CRT,M0SHA, 
§NPOINT,NOORPM, IASIGN, INPRPM,D,DN,P,EE,GG,EI,GAK,AM,ECC,AID,AIRO,BE 
§TA, GAMMA , BKMX, BKMY, BCMX , BCMY , XKMM , YKMM , XCMM , YCMM , BM , B I , QKXX , QKXY , Q 
5 KYY , QKYX , QCXX , QCXY , QCYY , QCYX , XXMK , XYMK , YYMK , YXMK , XXMC , XYMC , YYMC , YX 
&MC , KK , FDDF I X , BBB * BDB , BEB , BHB , BKB , BNB , BROB , QK , QC , QKP , QCP , QKF , QCF , 

$QKPF,QCPF, XKF , XCF , XKFF , XCFF , QKHD , QCHD , QKHDG , QCHDF , CT1,CT,CT2, 
§MT,MT1,MT2,AT,BT,DU,ET,HT,FT,GT, AA,BA,DA,EA.HA,FA,GA, GX,GY, 

§USV , USC , UBV , UBC , UTV , UTC , F 1 , ALFA, BCB , IPRINT, ID 

For data name description and units used, refer to Appendix B, Table XIV. 
Definition of Fortran Variables in Common Block. All namelist names ex- 
cept ID which is for input card identification only, may be found in the 
common block Fortran variable definition. 
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COMPUTER PROGRAM CAPABILITY AS 
CONTROLLED BY INPUT VARIABLES 

The values of certain input variables determine which area of computer capabil- 
ity will be activated or bypassed. The following is a description of the com- 
puter optional capability and related control varialbes. 

Rotor Mechanical Hysteresis Effects 

The hysteresis effects include 

1. In-phase transverse damping hysteresis forces and moments 

2. Out-of-phase transverse driving and damping hysteresis forces and 
moments 

3. In-phase damping torque 

There are six hysteresis coefficients which lead to the above hysteresis 
effects. These coefficients are: 

1. Transverse shear viscous hysteresis coefficient 

2. Transverse shear Coulomb friction hysteresis coefficient 

3. Transverse bending viscous hysteresis coefficient 

4. Transverse bending Coulomb friction hysteresis coefficient 

5. Torsional shear viscous hysteresis coefficient 

6. Torsional shear Coulomb hysteresis coefficient. 

Only by including any of the above coefficients for certain rotor sections (or 
all rotor sections) in namelist/OPTION/, will the related hysteresis portion of 
the computer program be activated. By inputing a hysteresis coefficient (or 
hysteresis coefficients), the hysteresis load as well as the dissipative hys- 
teresis torque will be computed and their effects on the rotor dynamic perform- 
ance included. 

Torsional Flexibility and its Control 

This parameter is used to simulate the torsional dynamic performance of a rotor. 
The incorporation of rotor torsional flexibility also makes it possible to in- 
clude the torsional hysteresis effects of the rotor. 

The optimum time step size used in a stable computation varies according to the 
absolute magnitude of local mass acceleration. The stable time step size is 
dominated by the maximum local acceleration in the rotor bearing system. To 
achieve a low computer time to real time ratio, it is desirable to have reason- 
able uniform acceleration rates among all local mass-load components. A large 
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local peak acceleration magnitude which may be caused by a very stiff, torsional 
rotor section among other sections results in a time consuming slow computation. 
Converting this very stiff section into a rigid section will in general lead to 
a computer time saving without affecting much of the computation accuracy. For 
this reason a rigid torsional section control parameter was incorporated in the 
program. With this parameter, the user can eliminate the torsional elasticity 
of certain rotor sections whose configurations are substantially stiffer than 
those of the others. An integer variable RIG (J) , is used to assign an artifi- 
tial rigidity of rotor section J by letting RIG (J) =1. If RIG (J) is not 
input, the built-in values of RIG (J) = 0, which considers actual torsional 
flexibility, will be used. 

Rotor Transverse Effects Due to Rotor Torsion 

For rotor sections having slope deflection, torsional loading will result in 
rotor transverse motion. For computations where this effect is desired, set 
IMT = 1 which overrides the builtin value of IMT = 0. 

Rotor Transverse Effects Due To Axial Loading 


The current rotor model does not have axial elasticity parameters, and the only 
dynamic effect of the axial loading is the rotor transverse motion. When 
IPP = 1 the rotor transverse motion effect due to axial loading is included. 
When the IPP value is not read in, the built-in value of IPP = 0 will be used 
to instruct the program to bypass the effects. 

Bearing In-Phase and Out-of-Phase Anisotropic 
Stiffness and Damping Force and Moment Coefficients 

In a specific rotor dynamics analysis, the effects of the coefficients can be 
included or deleted by using the appropriate input or the built-in default 
values of the coefficients. 

Mount In-Phase Anisotropic Stiffness and 
Damping Force and Moment Coefficients 

Similar treatment in applying or deleting the effects of the coefficients as 
stated in the bearing coefficient above may be used. 

Bearing Mass and Transverse Mass 
Moment of Inertia 


Choice of including or deleting the bearing mass and inertia effects may be 
exercised for the same reason discussed above for torsional flexibility. For 
a comparatively small bearing mass (or inertia) to bearing load ratio, use 
BM (K) = 0, or BI (K) = 0 (built-in value). This will instruct the program 
to bypass the bearing acceleration computation while maintaining reasonable 
computation accuracy. When a substantial BK (K) or BI (K) to bearing load 
ratio exists, use their actual values. 
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Additional Control Parameters 

There are other input control parameters such 

MET = 1 For using international units input and output 

MET = 0 For using English units input and output 

A complete list of input control parameters may be found in Section IV Final 
Verification of IBM 360/370 computer program. 

Other capability such as automatic restart, when the assigned integration time 
step (DT) is too large that the initial computation results in journal displace 
ments exceeding the bearing clearance. The computation is limited at present 
to a total of five restarts. In each restart, DT is reduced to one- fourth its 
previous value. 

COMPUTER PROGRAM OUTPUT AND ITS CONTROL 

The computer program output is in the form of printed output and graphic (CRT) 
plots as discussed below. The units of the input-output data may be selected 
by inputting MET = 1 for international units and MET = 0 for English units. 

Printout 


The printed output includes the following: 

1. Input Data Write-Out . Except that which is input through punched 
cards . 

2. Input Rotor Mass Data . Includes rotor local masses, transverse and 
polar mass moments of inertia, total rotor mass and polar mass moment 
of inertia, and location of the rotor mass center. 

3. Rotor Dynamic Startup Configuration . Includes deflections and loads 
for all rotor and bearing components. 


4. Computation Results Write-Out . The printout time interval can be 
approximately specified by the product of the values DP and IPRINT in 
seconds. 

5. Graphic (CRT) Output . This is an optional output controlled by CRT = 
1, or 0, and partially controlled by MOSHA = 1, or 0. When CRT = 1 
graphic output will be produced; CRT = 0 suppresses graphic output. 
The only type of graphic output control by MOSHA, in addition to CRT 
control variable, is the rotor mode- shape plot. When MOSHA = 1 and 
CRT =1 the rotor mode shape will be plotted, otherwise there will be 
no mode shape graphic output. A total of eight types of CRT graphs 



are provided where each presents a pictorial summary of certain rotor 
dynamic performance in supplementing the printed output. The time 
interval between successive points on the graph will be equal to the 
value of DP in seconds. The types of the graphic output are: 

a. Rotor mode shape graphs, one at or near each of the INPRPM (I). 

The phase angles of rotor deflection vectors are labelled at each 
rotor station. 

b. Rotor spin speed at station IASIGN versus time with the maximum 
time interval per graph equal to the product of DP and NPOINT. 

c. Rotor displacement whirl to spin velocity ratio versus time at 
rotor station IASIGN. NPOINT number of points will be included 
in each graph. 

d. Bearing force versus rotor spin speed for each of the support 
bearings. NPOINT number of points will be included in each graph. 

e. Bearing displacement versus rotor spin speed for each of the sup- 
port bearings. NPOINT number of points will be included in each 
graph . 

f. Maximum rotor deflection versus rotor spin speed, with the corres- 
ponding rotor station number labeled in the graph. NPOINT number 
of points will be included in each graph. 

g. Rotor deflection at station IASIGN versus rotor spin speed. 

NPOINT number of points will be included in each graph. 

h. Rotor orbital path for rotor station IASIGN. NPOINT number of 
points will be included in each graph 

The incorporation of graphic output capability in the computer program 
was provided. This output operation was not verified due to the absence 
of the specific graphic output facilities at Rocketdyne. 
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NONLINEAR BEARING STIFFNESS INPUT COEFFICIENTS DETERMINATION 


The nonlinear stiffness coefficients are provided to simulate a known or pre- 
dicted bearing stiffness characteristic. The steps in determining their 
coefficients are as follows: 

1. Load versus bearing displacement vector characteristic at a journal 
spin speed value ^ may be represented as in Fig. 57. 



Figure 57. A Bearing Force Versus Displacement Characteristics 



2. The bearing- load displacement characteristic for bearing stiffness, 
section K is now transformed into a corresponding chord stiffness as 
indicated by the dashed line (cd) in Fig. 58, 

BEARING 



Figure 58 . Bearing Stiffness vs Displacement Characteristics 

where (in Fig. 58) K is to be curve fitted by the formulation below, which 
is a displacement function of the nonlinear stiffness characteristics: 

K BiK L C BiK (V X Bi tY Bi ' »BiK ) ^BiK ( >/ X Bi* Y Bi ' P Bi K ) + E BiK. 

To determine the remaining constants, n d-i/ and „ in the nonlinear 
bearing formulation, an overall bearing chord stiffness curve ab 
for <f>i = ^Ki must be provided either from analytical prediction or 
experimental results. With this stiffness addition ac and bd for the 
spin speed <j)^ = ^Ki* the speed sensitive constant and can 

now be determined as follows 


N BiK 

^Ki " 

<f> .) 

oi _ ac 




K BiK 

CO 

r . 



^Ki 

- *oi> 

L N BiK + B BiK \\ 

/ X Bi * V Bi - w) 

bd 



SiK 


de 
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where 



Thus a combination of speed and displacement sensitive nonlinear bearing stiff- 
ness characteristic may be represented in various stiffness sections of a 
support bearing. 

PROGRAM SIZE CAPACITY 

The current computer program dimension size is compatible with the following 
maximums 


Rotor Stations (NS) 

\ 

Bearing Stations (NB) 

Nonlinear Bearing Stiffness sections (KK(K) ) = 3 

The current object program size including all necessary library and auxiliary 
requirements is approximately 199 K bytes. When the enlargement of the com- 
puter program capacity is required, the rule in adjusting the size of the 
dimension statement as indicated in Table VIII may be used. 


= 15 
= 6 
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TABLE VIII - RULES FOR MODIFYING COMPUTER PROGRAM CAPACITY 



Dimension Size 

Subscripted Variable 
Name or Group 

Current Computer 
Program Capacity 
NS = 15 
NB = 6 
KK(K) = 3 

Enlarged Computer Program 
Capacity, NS, NB, KK (K) 

Pertaining to Rotor 
Stations 

QM (15) 

QM(NS) 

C(15 , 15) 

C (NS ,NS) 


Etc. 

Etc. 

Pertaining to Rotor 

DD (14) 

DD (NS- 1) 

Section 

EE (14) 

EE (NS- 1) 


USV (14) 

USV (NS- 1) 


RIG (14) 

RIG (NS- 1) 


Etc. 

Etc. 

Pertaining to 

QKXX(6) 

QKXX(NB) 

Bearing Station 

BCB (6,3) 

BCB(NB,KK(K) ) 


BR0B (6 , 4) 

BR0B (NB , KK (K) +1 ) 


Etc. 

Etc. 

AA. . , in HYSSTA 
ij 

AA(84. 84) 

AA (4 (NS+NB) , 4 (NS+NB) 

CC. in HYSSTA 
1 

CC (84) 

CC (4* (NS+NB) ) 

YN. in MAIN 4 

l 

YN(84) 

YN (4* (NS+NB)) 

YN in FUND 

YN(198) 

YN(NNS) 

BD in FUNL 

BD (198) 

BD(NNS) 


Etc. 

Etc . 

Dimensioned Variables 
in HYSRKA, RUNKUT, 

Data NN/198/ 

Data NN/NNS/ 

ADAMLT 

Y ( 198) 

Y (NNS) 


F (198,7) 

F(NNS,7) 


A(198,4) 

A (NNS, 4) 


YP (198) 

YP (NNS) 


Etc. 

Etc. 


*NNS = 10NS+8NB 
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APPENDIX B 

DEFINITION OF FORTRAN VARIABLES 


The input Fortran variables definitions including their default values are 
described in Appendix C (Table XXI) and the Fortran variable except those in 
function ISIMEQ and subroutine PL0TBR are described and their units defined. 
ISIMEQ is a library subroutine in IBM 360/370. Its function is to solve a 
set of linear algebraic equations by inputting the coefficients and constants 
of the equation. The variables contained in ISIMEQ are purely mathematical 
notations which are not directly related to rotor dynamics analysis program- 
ming. The meaning of the variable names in PL0TBR are the same as those in 
MAIN and HYSMET. 

Fortran variable definitions in C0MM0N are described in Table IX. The varia- 
bles in each of the subroutines not covered in the C0MM0N table are defined 
in Tables X through XX according to the sequence they appear in each of the 
subroutines. All successive Fortran variable tables will only contain those 
variables not defined in the preceding tables. 



TABLE IX - DESCRIPTION OF ALL FORTRAN VARIABLES LISTED 

IN COMMON BLOCKS 


Variable 

Definition 

Units 

NS2 

= 2*NS 


NS3 

= 3*NS 


NS4 

= 4*NS 


NS5 

= 5*NS 


NS6 

- 6*NS 


NS7 

= 7*NS 


NS8 

= 8*NS 


NS9 

= 9*NS 


NS10 

= 10*NS 


NSM1 

= NS-1 


NSP1 

= NS'+l 


NS2P1 

= NS*2+1 


NS4P1 

= NS*4+1 


IP (or IQ 

Computer results printout frequency control 


IPRINT 

IPRINT x DP = computer output printout 
real time interval 


NN 

8*NB+10*NS; MM spacer in "common" in "fund" 


IB1 

Station number for the first bearing location 


IBNB 

Station number for the last bearing location 


INT 

Computation cycle indicator and control 
INT = 0 first time "fund" is called 
INT = 2 during "fund" calling process 


ITIM 

Startup control for integration 

ITIM = +1 restart with forward integration 
ITIM = -1 restart with backward integration 
ITIM = 0 continuing integration 


IUSE 

Interaction process and indicator 

IUSE = 1 completion of an integration 
process 

IUSE = 0 during integration < 


G 

Gravitational constant G = 386.088 

in. /sec 

PI 

it = 3.14159265358979324 
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TABLE IX - (Continued) 


Variable 

Definition 

Units 

Q 

Axial length between rotor station I and 
the last bearing station 

cm (inches) 

S 

Axial length between rotor station I and 
the first bearing station 

cm (inches 

QLL 

Axial length between the first and last 
bearing station 

cm (inches) 

QML0V 

Negative reciprocal of QLL 

-1 -1, 
cm (in. ) 

KKSPA(K) 

Number of nonlinear bearing stiffness sec- 
tions for bearing station K, KKSPA(K) is 
used as a spacer for KK(K) in the program 
where KK(K) is not used, while KK is used 
for other function 

dimensionless 

JBI(I) 

Bearing number at rotor station I 

JBI(I) = 0 for rotor station where no 
bearing exists 

dimensionless 

F(I) 

Starting rotor angular position for rotor 
station I 

degrees 

FD0T(I) 

Starting rotor angular spin speed for 

rotor station I ; 

rpm 

SHK(J) 

Average shear strain to actual shear strain 
ratio for rotor section J 

dimensionless 

QM(I) 

Total rotor mass at rotor station I 

2 

(kg-sec )/cm 
f(lb-sec 2 )/in.j 

QID(I) 

Total rotor transverse mass moment of inertia 
at rotor station I i 

' 2 
kg-cm-sec 

2 

(lb-in. -sec ) 

QIR0CI) 

Total rotor polar mass moment of inertia at 
rotor station I 

kg-cm-sec 2 

2 

(lb-in. -sec ) 

QME(I) 

Product of QM(I) and ECC(I) at rotor 

2 

kg-sec 

(lb-sec 2 ) 

F0STIFF(K) 

Nonlinear bearing stiffness at bearing 
station K 

kg/cm (lb/in.) 

Z(I) 

Rotor axial length measured from rotor 
stations 1 to (I) 

cm (inches) 

QZ(I) 

Q less the amount Z(I) 

cm (inches) 

QK(I) 

Rotor-to-casing in-phase stiffness force 
coefficient at rotor station I 

Newtons/cm 

(lb/in.) 
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TABLE IX - (Concluded) 


Variable 

Definition 

Units 

QC ( I ) 

Rotor-to-casing in-phase damping force 
coefficient at rotor station I 

(Newton-sec)/cm 

(lb-sec)/in. 

YN (M) 

Startup rotor and bearing displacement and 
slope array M = 4* (NS + NB) 

inches, radians 

YNSPA(M) 

"Common” spacer for YN(M) 


INPRPM(N) 

Number of rotor spin-speed rpm one for each 
rotor mode shape CRT plot 


C(I,I1) 

Rotor linear deflection at station 11 from the 
straight line joining the first and last bear- 
ing centers due to unit transverse loading at 
station I 

in. /lb 

B(I > 11) 

Rotor linear deflection at station 11 from the 
straight line joining the first and last bear- 
ing centers due to unit transverse moment 
loading at station I 

in. /lb-in. 

TF ( 1 , 1 1 ) 

Rotor slope deflection at station 11 from the 
straight line joining the first and last bear- 
ing centers due to unit transverse loading at 
station I 

radians/lb 

TM(I , 11) 

Rotor slope deflection at station 11 from the 
straight line joining the first and last bear- 
ing centers due to unit transverse moment 
loading at station I 

radians/lb-in. 
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TABLE X - DEFINITION OF FORTRAN VARIABLES USED IN MAIN PROGRAM 


Variable 

Definition 

Units 

A 

180/tt conversion factor between angular 
degrees and radians 

degrees/radian 

H 

30/it conversion factor between rpm and 
radians/sec 

rev-sec 

r.adian-min 

I,J,M 

Subscripts used in dimensioned variable 


V 

0.5/tt conversion constant 


IC 

CRT variable accumulator subscript 


II 

CRT mode shape plot counter 


IP 

Printing interval counter 


IS 

Output printing interval counter 


KA 

Punch card ID sequence specification 


KB 

Punch card block ID sequence specification 


R0(I) 

Rotor displacement vector 

inches 

SP(N) 

Transfer of YNN(I) to punch 6 values per 
card 

inch, in. /sec 

TT(N1) 

Subscripted T for CRT data accumulation 

seconds 

XB(K) 

Bearing X-displacement at bearing station K 

inches 

XX (I) 

Rotor X-displacement at rotor station I 

inches 

YB(K) 

Bearing Y-displacement at bearing station K 

inches 

YY(I) 

Rotor Y-displacement at rotor station I 

inches 

DDA 

Printing interval control variable 

seconds 

INS, JNS,MNS 

Subscripts 


REV (I) 

Rotor spin revolution at rotor station I 
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TABLE X - ( Continued) 


Variable 

Definition 

Units 

RPM(I) 

Rotor spin speed at rotor station I 

rpm 

TSA 

Initial "T" saved for restart with 
10 percent DT 

seconds 

XBM(K) 

Bearing slope in XZ-plane at bearing 
station K 

radians 

XXT(N) 

Rotor X- displacement component printing 
time-step accumulation (IC) array at rotor 
station "IASIGN," for CRT use 

inches 

YBM(K) 

Bearing slope in YZ-plane at bearing 
station K 

radians 

YYT(N) 

Rotor Y-displacement component printing 
time-step accumulation (IC) array at rotor 
station "IASIGN," for CRT use 

inches 

BRGR(IC,K) 

Bearing displacement vector printing time- 
step accumulation (IC) array at bearing 
station K, for CRT use 

inches 

F0RC(IC,K) 

Bearing force vector printing 
time-step accumulation (IC) array at 
bearing station K, for CRT use 

pounds 

I ERR 

Indicator from subroutine HYSRKA 
IERR = 0 solution is valid 
IERR = 1 solution is invalid or 
integration time step 
reaches zero 


I10S,I2NS 
I3NS, I4Ns 
I5NS, I6NS 
I7NS.I8NS 
I9NS, J2NS 
J3NS, J4NS 
J5NS, J6NS 
J7NS, J9NS 

Subscripts used in dimensioned variable 


M0R0 ( I ) 

Mount displacement vector 

inches 

M0SQ 

Square of mount displacement vector 

. 2 
m. 

M2NS.M3NS, 

M5M1 

Subscripts used in dimensioned variables 
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TABLE X - (Continued) 


Variable 

Definition 

Units 

R0MM(I) 

Mount slope vector 

inches 

R0SQ 

Square of rotor displacement at a rotor 
station 

. 2 

in. 

RPMM(IC) 

Rotor spin speed array at rotor station 
IASIGN for A series of time step; A CRT 
plotting variable 

rpm 

SL0P(I) 

Rotor slope vector at rotor station I 

radians 

XM0M(I) 

Mount moment in XZ-plane at rotor station I 

lb-in. 

YM0M(I) 

Mount moment in YZ-plane at rotor station I 

lb- in. 

BRGR0 (K) 

Bearing displacement vector at bearing 
station K 

inches 

BSLR0(K) 

Bearing slope vector at bearing station K 

radians 

DTAVE 

Average integration step time 

seconds 

I10SB 

A subscript used in dimensioned variables 


M0F0R(K) 

Mount force vector at station K 

pounds 

PHAR0(I) 

Rotor displacement vector phase angle 

degrees 

ROMAX(IC) 

Rotor maximum displacement vector printing 
time-step accumulation (IC) array at rotor 
station "IASIGN" for CRT use 

inches 

R0STA(IC) 

Rotor displacement vector printing time-step 
accumulation (IC) array at rotor station 
"IASIGN" for CRT use 

inches 

TSAVE 

Time saved for computing average integration 

seconds 

WHIRR (I) 

Rotor displacement whirl frequency 

rpm 

XBDOT(K) 

Bearing displacement X-velocity at bearing 
station K 

in. /sec 

XBFOR(K) 

Bearing Y-force at bearing station K 

pounds 

XBMON(K) 

Bearing moment in XZ-plane at bearing station K 

lb- in. 








TABLE X - (Continued) 


Variable 

Definition 

Units 

XMF0R(K) 

Mount X- force at bearing K 

pounds 

YBD0T(K) 

Bearing displacement Y-velocity at bearing 
station K 

in. /sec 

YBF0R(K) 

Bearing Y-force at bearing station K 

pounds 

YBMOM(K) 

Bearing YZ-plane moment at bearing station K 

lb-in. 

YMF0R(K) 

Mount Y-force at bearing station I 

pounds 

BGPHAS(K) 

Bearing displacement vector phase angle at 
bearing station K 

degrees 

BRFOPH(K) 

Bearing force vector phase angle at bearing 
station K 

degrees 

BRGFOR(K) 

Bearing force vector at bearing station K 

pounds 

BSPHAS(K) 

Bearing slope vector phase angle at bearing 
station K 

degrees 

ISTATN(IC) 

Rotor station time-step array for maximum 
rotor displace at each time-step for CRT 
graph 


I10S2B, I10S3B 
I10S4B, I10S5B 
I10S6B, I10S7B 
14NSNB 
14NS2B 
14NS3B 

Subscripts used in dimensioned variables 


M0F0PH(K) 

Mount force vector phase angle at bearing 
station K 

degrees 

M0PHAS(K) 

Mount displacement vector-phase angle at 
bearing station K 

degrees 

M0WHIR(K) 

Mount displacement (bearing mass) whirl fre- 
quency also mount whirl to rotor spin speed 
(bearing mass) whirl frequency ratio at 
bearing station K 

radians/sec 

dimensionless 

PHAR0S ( I ) 

Rotor slope vector phase angle at rotor 
station I 

degrees 

PHASMM(K) 

Mount slope vector phase angle at bearing 
station K 

degrees 

SLOPSQ 

Rotor slope vector square at a rotor station 

2 

(radians) 

WHRATI(IC) 

Rotor displace whirl to rotor spin velocity 
ratio time-step array at rotor station 
"IASIGN" 

dimensionless 
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TABLE X - (Concluded) 


Variable 

Definition 

Units 

WHRAT0(I) 

Rotor displacement whirl-to-rotor spin velocity 
ratio at rotor station I 


WHRVL0(I) 

Rotor displacement whirl velocity at rotor 
station I 

radians/sec 

WHSL0P (I) 

Rotor slope whirl velocity at rotor station I 

rpm 

XBMD0T(K) 

Bearing slope velocity in XZ-plane at bearing 
station K 

radians/sec 

YBMD0T (K) 

Bearing slope velocity in YZ-plane at bearing 
station K 

radians/sec 

YNNSAV ( I ) 

Saved YNN(I) data from a continued fun when 
"C0NTIN=1" 


IR 

Number of restart control variables 


rH 

Subscripts used in dimensioned variables 


ACA 

Journal displacement vector at a bearing 
station 

inches 
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TABLE XI - DEFINITION OF FORTRAN VARIABLES USED IN SUBROUTINE 

HYSREA 


Variable 

Definition 

Units 

I,J,K 

Subscripts used in dimensioned variables 


API, AP2 , 
AP3, AP4 

Variables for detecting the axial load unbalance 
with dimension corresponding to that of AA(I), 
BA(I) , DA (I) and EA(I) 


NB4 

Four times number of bearing 


MUST 

Name of one of the two name lists used to read 
input data. Variables contained in this namelist 
must be input to the program at all times 


0PTI0N 

Name of one of the two namelists used to read 
input data. The variables contained in this 
namelist may be read in total or in part as 
desired. 
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TABLE XII - DEFINITION OF FORTRAN VARIABLES USED IN SUBROUTINE 

HYSWRI 


Variable 

Definition 

Units 

I,J,K 

Subscripts used in dimensioned variables 


U 

it divided by 30. A conversion constant 

radians-min 


between rotational speed in rpm and that in 
radians/sec 

rev-sec 

V 

ir divided by 180. A conversion constant 
between angular displacement in degrees and 
radians 

radians/degree 

K1 

Maximum value of the subscript KI for 
BR0B(K,KI) for each of the nonlinear 
stiffness bearings 












TABLE XIII - DEFINITION OF FORTRAN VARIABLES USED IN SUBROUTINE 

HYSWME 


Variable 

Definition 

Units 


Subscripts used in dimensioned variables 


U,V 

Same as those defined in Table XII 


AF 

Conversion constant from pounds to Newtons 

Newtons/lb 

IC,II 

Same as those defined in Table X 


KI 

KI=KK(K) 

Number of stiffness sections for each of the 
nonlinear stiffness bearings 


K1,K2 

Maximum value of the subscript KI for 
BROB(K,KI) for each of the nonlinear 
stiffness bearings 


ADN 

Material weight density conversion constant 
from lb/in.^ to kg/cm^ 


AIN 

Equivalent of cm to inch 


AFIN 

Conversion constant from lb-in. to Newton-cm 


AFIN 2 

2 

Conversion constant from lb-in. to 
Newton- cm^ 


AF0IN 

Conversion constant from lb/in. to Newtons/cm 


AINER 

2 2 

Conversion constant from lb-in. to kg-cm 


AMASS 

Conversion constant from pounds to kg 


AF0IN2 

Conversion constant from lb/in. to 
Newtons/cm 

i 

■ 
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TABLE XIV - DEFINITION OF FORTRAN VARIABLES USED IN SUBROUTINE 

HYSMET 


Variable 

Definition 

Units 

A,H 

Same as that defined in Table X 


I.J 

Subscripts used in dimensioned variables 


V 

Same as that defined in Table X 


AIN 

Same as that defined in Table XIII 


REV (I) 

Same as that defined in Table X 


XXT(IC) 

Rotor X-displacement component printing time- 
step accumulation (IC) array at rotor station 
"IASIGN," for CRT use 

cm 

YYT(IC) 

Rotor Y-displacement cpmponent printing time- 
step accumulation (IC) array at rotor station 
"IASIGN,” for CRT use 

cm 

F0RC(IC,K) 

Bearing force vector printing time-step 
accumulation (IC) array at bearing station K, 
for CRT use 

Newton 

BRGR(IC,K) 

Bearing displacement vector printing time-step 
accumulation (IC) array at bearing station K, 
for CRT use 

cm 

R0SQ,RPMM ' 

SL0P,SQRT 

IlOSB 


Same as that defined in Table X 


BRGR0 (K) 

Bearing displacement vector at bearing 
station K 

inches 

ROMAX(IC) 

Maximum rotor displacement vector (at a time 
point) printing time-step accumulation (IC) 
array for CRT use 

cm 

R0STA(IC) 

Rotor displacement vector printing time^step 
accumulation (IC) array at rotor station 
"IASIGN" for CRT use 

cm 
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TABLE XIV - (Concluded) 


Variable 

Definition 

Units 

XBD0T, XBF0R, 

XBM0M,XMF0R, 

XBD0T,YBF0R, 

YBM0M,YMF0R, 

BRGF0R , BSPHAS , 

ISTATN 

Same as those 

defined in Table X 


I10S2B,I10S3B, 
I10S4B, I10S5B, 
I10S6B,I10S7B, 




M0F0PH ,M0PHAS , 
M0WHI R , PHAR0S , 
PHASN1M,SL0PSQ 




WHRATI 

WHRVL0 

WHSL0P 

XBMD0T 

YBMD0T 
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TABLE XV - DEFINITION OF FORTRAN VARIABLES USED IN SUBROUTINE 

HYSINF 


Variable 

Definition 

Units 

I,J,K 

Subscripts used in dimensioned variables 


RA(J) 

Rotor outside diameter to inside diameter 
ratio for rotor section (J) 

dimensionless 

SZ(I) 

Z-coordinate (axial) value from Ith rotor 
station to first bearing station 

inches 

ZQ Cl) 

Z-coordinate of rotor station I less that at 
the last bearing station 

inches 

F0F(I,I1) 

Reactional transverse rotor force at station 
11 due to unit force application at rotor 
station I 

pounds 

F0M(I,I1) 

Reactional transverse rotor force at rotor 
station 11 due to unit moment application at 
rotor station I 

lb-in. 

IB2 

Last bearing rotor station number 


M0F (1,11) 

Reactional rotor transverse moment at rotor 
station I due to unit force application at 
rotor station I 

pounds 

M0M(I,I1) 

Reactional rotor transverse moment at rotor 
station 11 due to unit moment application at 
rotor station I 

lb-in. 

R0F(I,J) 

Rotor transverse deflection at rotor station J 
relative to the straight line joining first 
and last bearing centers due to load at rotor 
station I 

in. /lb 

R0M(I,J) - 

Rotor transverse deflection at rotor station J 
relative to the straight line joining first 
and last bearing centers due to moment at 
rotor station I 

in. /(lb-in.) 

OLE I (J) 

'■ c . . 

Rotor sectional length divided by the product 
of Young's modulus of elasticity and sectional 
area moment of inertia for rotor section J 

1/lb-in.) 

QZ0L( 1) 

Negative of ZA(I) divided by the span between 
first and last bearing 

dimensionless 
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TABLE XV 


(Concluded) 


Variable 

. 1 

Definition 

Units 

SZ0L(I) 

SZ(I) divided by the span between first and 
last bearing 

dimensionless 

ZQ0L(I) 

ZA(I) divided by the span between first and 
last bearing 

dimensionless 

AFALEI (I) 

Sum of SHERGA(J) and cube of rotor sectional 
length divided by three times EI(J) for rotor 
section J 

in. / lb 

SHERGA(J) 

Rotor sectional length divided by the product 
of modulus of shear rigidity, cross-sectional 
area and inverse of shear strain factor (GAK(J)) 
for rotor section J 

in. /lb 

SQLZEI(J) 

Rotor sectional length square divided by 2EI (J) 
for rotor section J 

lb' 1 








TABLE XVI - DEFINITION OF FORTRAN VARIABLES USED IN SUBROUTINE HYSSTA 


Variable 

Definition 

Units 

I , J,K,M 

Subscripts 


U 

[ 

4/3, a constant 

dimensionless 

V 

180 divided by tt, radian to angular degrees 
conversion constant 

degrees/radian 

W 

tt divided by the product of 128* gravitational 
constant 

' - 

2 

sec /m. 

AA(I ,11) 

Matrix coefficients for solution of startup 
rotor deflection configuration 

dimensionless 

CC (I ) 

Matrix coefficient constant terms for solution 
of startup rotor deflection configuration 

inches , radians 

CF (I) 

Square of initial rotor spin speed times rotor 
mass at rotor station I 

lb/in. 

CG 

Center of rotor mass from rotor station 1 

inches, cm 

CM 

Square of initial rotor spin speed times the 
difference of rotor transverse mass moment of 
inertia less the polar mass inertia at a rotor 
station 

lb-in. 

FC 

Product of rotor mass eccentricity vector, 
square of initial rotor spin speed and rotor 
mass at a rotor station 

pounds 

FF 

A variable whose value is used as a scale 
factor for determinant evaluation in the 
"ISIMEQ" library function FF=0 is used here 


FG 

Sum of initial rotor spin angular position 
and rotor eccentricity phase angle at a rotor 
station 

radians 

FI 

Initial rotor spin angular position 

radians 

IA 

A variable having one-dimensional erasable 
array at least equal to the number of rows of 
matrix AA(1, 11) . for -library function "ISIMEQ” 
use 
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TABLE XVI - (Continued) 




Variable 

Definition 

Units 

IK, JK, JS 

Subscripts used in dimensioned variables 


KL 

Integer variable name for the library function 
"ISIMEQ". In return to the calling program KL 
serves as an indicator of the status of solu- 
tion, i.e., KL=1 successful solution 

KL=3 singular matrix detected 


MC 

Product of square of initial rotor spin speed, 
value of rotor transverse mass moment of 
inertia less polar mass moment of inertia and 
the misalignment angle of the rotor mass moment 
of inertia 

lb-in. 

QW Cl) 

Rotor mass at rotor station I 

pounds , kg 

R0(I) 

Rotor displacement vector at rotor station I 

inches , cm 

XB(K) 

Bearing X-displacement at bearing station K 

inches 

XM(K) 

Mount X-displacement at bearing station. K 


XS(I) 

Rotor slope in XZ-plane at rotor station I 

radians 

XX (I) ' \ 

Rotor X-displacement at rotor station I 

inches 

YB (K) 

''Bearing Y-displacement at bearing station K 

inches 

YM(K) 

Mount Y-displacement at bearing station K 

inches 

YN (N) 

Computed rotor and bearing startup displace- 
ments and slopes 

inches , radians 

YS (I) 

Rotor slope in YZ-plane at rotor station I 

radians 

YY(I) 

Rotor Y-displacement at rotor station I 


C0N (N, 1) 

Computed rotor and bearing startup displace- 
ments and slopes 

inches, radians 

DDL (J) 

Sum of square of rotor outside diameter, square 
of rotor inside diameter and four thirds 
square of rotor sectional length for rotor 
station J - - 

. 2 
in. 
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TABLE XVI - (Continued) 


Variable 

Definition 

Units 

DD2 (J) 

Square of rotor outside diameter for rotor 
section (J) 

. 2 
in. 

FAA 

Sum of the rotor initial starting spin angular 
displacement and rotor mass eccentricity phase 
angle at a rotor station 

radians 

FXC(I) 

Negative of cosine of eccentricity vector 
angular displacement times "FC" at rotor 
station I 

pounds 

FXX(I) 

X-force gradient for X-displacement due to 
mass and related nonbearing stiffness and 
damping coefficients at rotor station I 

lb/in. 

FXY(I) 

X-force gradient for Y-displacement due to 
related nonbearing stiffness and damping 
coefficients at rotor station I 

lb/in. 

FYC(I) 

' ' ' ' 1 

Negative of sine of eccentricity vector 
angular displacement times "FC" at rotor 
station I 

pounds 

INB, JKS, JNB, 
JNS, J2S, J3S, 
MNS . 

Subscripts used in dimensioned variables 


MXC(I) 

Negative cosine of misalignment vector 
angular displacement times MC 

lb-in. 

MXX(I) 

XZ -plane moment gradient for XZ-plane slope 
due to mass inertia and related nonbearing 
stiffness and damping coefficients at rotor 
station I 

lb-in. /radian 

MXY(I) 

XZ-plane moment gradient for YZ-plane slope 
due to related nonbearing stiffness and 
damping coefficients at rotor station I 

lb- in. /radian 

MYC(I) 

Negative; sine of misalignment vector angular 
displacement times MC 

lb-in. 

QKB(K) 

Startup nonlinear bearing stiffness at bearing 
station K 

1 =: 

lb/in. 

i - 
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TABLE XVI - (Continued) 


Variable 

Definition 

Units 

QL2 (J) 

Square of rotor sectional length for rotor 
section J 

. 2 

in . 

SIN (A) 

Sine function with argument "A" in radians 


SUM 

Total rotor mass moment about rotor station I 

lb-in. 

SBM(K) 

Bearing slope at bearing station K 

radians 

XMM(K) 

Mount slope at bearing station K 

radians 

XXC (K) 

Average in-phase bearing damping force 
coefficient of the corresponding nonisotropic 
X- and Y-components at bearing station K for 
startup computation purpose 

lb-sec/in. 

XXK(K) 

Average in-phase bearing stiffness force 
coefficient of the corresponding nonisotropic 
X- and Y-components at bearing station K for 
startup computation purpose 

lb/in. 

XYC(K) 

Average out-of-phase bearing damping force 
coefficient of the corresponding nonisotropic 
X- and Y-components at bearing station K for 
startup computation purpose 

lb-sec/in . 

XYK(K) 

Average out-of-phase bearing stiffness for 
coefficient of the corresponding nonisotropic 
X- and Y-components at bearing station K for 
startup computation purpose 

lb/in. 

YBM(K) 

Bearing slope in YZ-plane at bearing station K 

radians 

AFIN 

A conversion constant from lb-in. to Newton-cm 

Newton-cm/ 

lb-in. 

BCMM(K) 

Mount damping moment coefficient at bearing 
station K 

lb-in. -sec/ 
radian 

BKMM(K) 

Mount stiffness moment coefficient at bearing 
station K 

lb-in. /radian 
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TABLE XVI - (Continued) 


Variable 

Definition 

Units 

I2NB.I2NS, 

I3NB.I3NS, 

I4NS.JK2S 

JK3S,J2NB, 

J2NS,J3NB, 

J3NS,M2NS, 

M3NS 

Subscripts used in dimensioned variables 


NS41 

One plus four times number of rotor stations, 
used as a D0 L00P index 


QIDW(I) 

Rotor transverse mass moment of inertia at 
rotor station I 

.. . 2 
lb-in. 

kg -cm^ 

R0MM(K) 

Mount slope vector at bearing station K 

radians 

R0SL(I) 

Rotor slope vector at rotor station I 

radians 

WEIT 

Total rotor mass 

pounds, kg ! 

XBD0(K) 

Bearing linear X- velocity at bearing station K 

in. /sec 

XMD0(K) 

Mount linear X-velocity at bearing station K 

in. /sec 

XXCM 

Average in-phase bearing damping moment 
coefficient of the corresponding nonisotropic 
X- and Y-components at bearing station K for 
startup computation purpose 

(lb-in. -sec)/ 
radian 

XXKM 

Average in-phase bearing stiffness moment 
coefficient of the corresponding nonisotropic 
X- and Y-components at bearing station K for 
startup computation purpose 

(lb-in.)/ 

radian 

XYCM 

Average out-of-phase bearing damping moment 
coefficient of the corresponding nonisotropic 
X- and Y-components at bearing station K for 
startup computation purpose 

(lb-in. -sec)/ 
radian 

XYKM 

Average out-of -phase bearing stiffness moment 
coefficient of the corresponding nonisotropic 
X- and Y-components at the bearing station K 
for startup computation purpose 

(lb-in. ) /sec 
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TABLE XVI - (Continued) 


P 


Variable 

Definition 

Units 

YBD0 

Bearing linear Y-velocity at bearing station K 

in. /sec 

YMD0 

Mount linear Y-velocity at bearing station K 

in. /sec 

ZS0L 

The value of rotor axial (Z) coordinate less 
that of the first bearing then divided by the 
span between first and last bearing 

dimensionless 

AMASS 

Conversion constant from pounds to kg 

kg/ lb 

AMIN2 

2 2 

Conversion constant from lb-in. to kg-cm 

kg-cm 2 /lb-in. 2 

ATAN 2 (Y , X) 

Arc tangent function using Y, and X argument 

radians 

BRGR0(K) 

Bearing displacement vector at bearing 
station K 

inches , cm 

BSLR0(K) 

Bearing slope vector at bearing station K 

radians I 

COSFA 

Cosine of rotor mass eccentricity vector 
displacement angle at a rotor station 

dimensionless 

COSFG 

Cosine of rotor mass inertia misalignment 
displacement angle at a rotor station 

dimensionless 

DDPLD (J) 

Sum of the square of rotor outside diameters 
and square of rotor inside diameter for rotor 
section J 

. 2 
in. 

IBINS,I4NSB 

Subscripts used for dimensioned variables 


NS4NB 

Total number of bearings plus four times 
total number of rotor stations 

dimensionless 

QIR0W (I) 

Rotor polar mass moment of inertia at rotor 
station I 

lb-in 2 
kg- cm^ 

QMASS 

Total rotor mass 

2 

lb-sec /in. 

SINFA 

Sine of rotor mass eccentricity vector displace- 
ment angle at a rotor station 

dimensionless 

SINFG 

Sine of rotor mass inertia misalignment dis- 
placement angle at a rotor station 
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TABLE XVI - (Continued) 


Variable 

Definition 

Units 

XBF0R(K) 

Bearing X-force at bearing station K 

pounds .Newtons 

XBIM0(K) 

Bearing mass moment of inertia XZ-plane 
moment at bearing station K 

lb-in. , 
Newton-cm 

XBMD0(K) 

. 

Bearing slope XZ-plane velocity at bearing 
station K 

radians/sec 

XBMF0 (K) 

Bearing mass inertia X-force at bearing 
station K 

pounds, Newtons 

XBM0M(K) 

XZ-plane bearing moment at bearing station K 

lb-in. , 
Newton-cm 

XMF0R(K) 

Mount X-force at bearing station K 

p oun d s , N e w t on s 

XMMD0(K) 

Mount XZ-plane slope velocity at bearing 
station K 

radians/sec 

XMMOM(K) 

Mount XZ-plane moment at bearing station K 

lb-in. , 
Newton-cm 

YBF0R(K) 

Bearing Y-force at bearing station K 

pounds , Newtons 

YBIM0 (K) 

Bearing mass moment of inertia YZ-plane 
moment at bearing station K 

lb- in., 
Newton-cm 

YBMD0 (K) 

Bearing slope ZY-plane velocity at bearing 
station K 

radian/sec 

YBMF0(K) 

Bearing mass inertia Y-force at bearing 
station K 

pounds , Newtons 

YBM0M(K) 

YZ-plane bearing moment at bearing station K 

lb-in., 

Newton-cm 

YBF0R(K) 

Bearing Y-force at bearing station K 

pounds, Newtons 

YMMD0(K) 

Mount YZ-plane slope velocity at bearing 
staion K 

radians/sec 

YMM0M(K) 

Mount YZ-plane moment at bearing station K 

lb-in. , 
Newton-cm 

BRPHAS(K) 

Bearing displacement vector phase angle at 
bearing station K 

degrees 
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TABLE XVI - (Continued) 


Variable 

Definition 

Units 

BSPHAS(K) 

Bearing slope vector phase angle at bearing 
station K 

degrees 

FD0TSQ 

Square of starting rotor spin speed 

radians/sec 

IBNBNS 

A subscript used in dimensioned variables 


ISIMEQ 

Library function name for solution of 
simultaneous equations 


I4NS2B, I4NS3B 

Subscripts used in dimensioned variables 


MOPHAS (K) 

Mount displacement vector phase angle at 
bearing station K 

degrees 

M0UNR0 (K) 

Mount displacement vector at bearing station K 

inches , cm 

NS42NB 

Sum of four times total number of rotor 
stations and two times total number of bearing 
stations 


NS43NB 

Sum of four times total number of rotor 
stations and three times total number of 
bearing stations 


NS44NB 

V 

Sum of four times total number of rotor 
stations and four times total number of 
bearing stations 


PHAR00 ( I ) 

Rotor displacement vector phase angle at rotor 
station I 

degrees 

PHAR0S (I) 

Rotor slope vector phase angle at rotor 
station I 

degrees 

PHASMM(K) 

Mount slope vector phase angle at bearing 
station K 

degrees 

P0LARA 

Total rotor polar mass moment of inertia 

i ' ~ 

lb-in. ^ ^ 
lb-in. sec 
kg-cm-sec^ 
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TABLE XVI - (Concluded) 


Variable 

Definition 

Units 

PRIMAS 

Printout of rotor mass properties control 
variable 

PRIMAS=0, delete printing 
PRIMAS=1, print the properties 


PRISTA 

Print of startup rotor dynamic deflection 
control variable 

PRISTA=0, delete printing 
PRISTA=1, print the deflection 


QLDNDD 

4/tt times the rotor section mass for a rotor 
section 

pounds 

QILDND(J) 

Rotor sectional mass divided by 2 * 386.088 
for rotor section J 

lb-sec^/in. 

Q6LDND (J) 

i 

Rotor sectional mass divided by 32 * 386.088 
for rotor section J 

lb-sec /in. 



TABLE XVII - DEFINITION OF FORTRAN VARIABLES USED IN SUBROUTINE HYSRKA 


Variable 

Definition 

Units 

A(NN,I1), 

Integration processing variables 

First Subscript = number of input and output 

variables to be integrated 

Second Subscript = number of internal 

processing steps 


H 

Integration time step could have time or other 
units 


I,J 

Subscripts used for dimensioned variables 


K 

Indicator for overflow conditions 


L 

A flow procedure control constant 


N 

Total number of variables to be integrated 


NN 

Maximum allowable number of variables to be 
integrated 


NS 

A flow control variable 


XB 

Substitute time variable 


YB 

Substitute variable for Y(NN) 


NERR 

Integration validity indicator 
NERR = 0 solution is valid 
NERR = 1 solution is invalid either due to 
N is invalid or H has gone to zero 


ISTFLG 

A flow process control variable 


OVERFL 

A library subroutine which provides for testing 
for an exponent overflow or underflow in real 
(floating point) operations 



142 







TABLE XVIII - DEFINITION OF FORTRAN VARIABLES USED IN SUBROUTINE RUNKUT 


Variable 

Definition 

Units 

A(NN, 13) 

Derivatives supplied by "FUND" for integrating 
variables. Their units vary according to the 
units used for dependent and independent 
variables. 


H 

Time step or step of other type of independent 
variable used 


I 

A subscript for a dimensioned variable 


N 

Number of variables to be integrated 


V(I) 

An integration variable array 


X 

An integration time or other independent 
variables 


Y 

Variables to be. integrated 


XI 

Alternate time (or other variable) step variable 


X2, X3 

Alternate time (or other variable) in integration 


YY(NN) 

Solution of integration variable at end of time 
(or other variable) step 









TABLE XIX - DEFINITION OF FORTRAN VARIABLES USED IN SUBROUTINE ADAMLT 


Variable 

Definition 

Units 

A (NN ,11) 
F (NN, I 2) 

An integrating variable array in the process 
of computation, 

NN = number of integration variable 

11 = number of internal processing levels 


H 

Time (or other variable) step 


I,K 

Subscripts used in dimensioned variables 


N 

Number of integration variables 


S 

A flow control variable 


T, U, V, W 

Equivalent variable names for testing inter- 
mediate integration variables 


X 

Time or other independent variables in 
integration 


Y 

Integration variables 


HH 

One- twenty- fourth of time (or other variables) 
step 


IC 

A flow process control variable 


YC (NN, 13) 

Corrector integration solution variables 


YP (NN ,14) 

Predictor integration solution variables 


NERR 

Same as that defined in Table 


ISTFLG 

A flow control flag 


OVERFL 

Same as that defined in Table 
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TABLE XX - DEFINITION OF FORTRAN VARIABLES USED IN SUBROUTINE FUND 


Variable 

Definition 

Units 

I » J»K 

Subscripts used in dimensioned variables 


AR(I) 

Equivalent polar mass moment of inertia at 
rotor station I 

lb-in. ^ 

EX 

Rotor XZ-plane slope differential for a rotor 
section 

radians 

EY 

Rotor YZ-p lane slope differential for a rotor 
section 


FG 

Rotor mass inertia misalignment angular dis- 
placement angle at a rotor station 

radians 

FX(I) 

Combined rotor transverse elastic and 
hysteresis X-force at rotor station I 

pounds 

FY (I) 

Combined rotor transverse elastic and 
hysteresis Y-force at rotor station I 

pounds 

HX 

Product of rotor sectional length and its 
XZ-plane slope less XMX for a rotor section 

inches 

HY 

Product of rotor section length and its ,YZ- 
plane slope less YMY for a rotor section 

inches 

IJ 

Subscript for a dimensioned variable 


IQ 

Dummy spacer in common 


IR 

A variable in torsional rigid section control 


JJ,Ki 

A subscript used in dimensioned variables 


K4 

Number of stiffness sections for a nonlinear 
stiffness bearing 


MX (I) 

Combined rotor transverse elastic and 
hysteresis XZ-plane moment at rotor station I 

lb-in. 

MY (I) 

Combined rotor transverse elastic and 
. hysteresis YZ-plane jnoment at rotor station I 

lb-in. 

PP(I) 

Total axes loading at rotor station I 

pounds 

UE 

Viscous rotor bending hysteresis coefficient 
divided by rotor Young's modulus of elasticity, 
for a rotor section 

seconds 
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TABLE XX - (Continued) 


Variable 

Definition 

Units 

UG 

Coulomb friction rotor shear hysteresis 
coefficient divided by rotor shear modulus 
of rigidity at a rotor section 

dimensionless 

XB(K) 

Bearing X-displacement at bearing station K 

inches 

YB(K) 

Bearing Y-displacement at bearing station K 


YN (NN) 

Incoming variables into "fund" from which the 
derivatives are generated 

inches, in. /sec, 
radians , 
radians/sec 

ACA 

Bearing displacement vector at a bearing 
station 

inches 

BXL(I) ,BXR(I) 

Combine rotor transverse XZ-plane, viscous 
and Coulomb friction hysteresis bending 
moment at the immediate left and right of 
rotor station I, respectively 

lb-in. 

BYL(I) , BYR(I) 

Combined rotor transverse YZ-plane, viscous 
and Coulomb friction hysteresis bending 
moment at the immediate left and right of 
rotor station I, respectively 

lb-in. 

FAA 

Rotor mass eccentricity vector angular 
displacement 

radians 

FDD (I) 

Rotor spin acceleration at rotor station I 

radians/sec^ 

FDL,FDR 

Absolute value of the rotor transverse force 
whirl to rotor spin velocity ratio at the 
immediate left and right of a rotor station, 
respectively 

dimensionless 

FHX(I) , FHY ( I ) 

Combined rotor transverse viscous and Coulomb 
friction hysteresis X- and Y-force at rotor 
station I, respectively 

pounds 

FLS,FRS 

Rotor transverse elastic force vector at the 
immediate left and right of a rotor station, 
respectively 

pounds 


146 








TABLE XX - (Continued) 


Variable 

Definition 

Units 

FWL,FWR 

Rotor force spin-whirl velocity difference to 
its absolute value ratio at the immediate left 
and right of a rotor station, respectively 

dimensionless 

FXL(I) ,FXR(I) 

Rotor elastic X-force to the immediate left 
and right of a rotor station I, respectively 

pounds 

FXX 

An X-force function 


FYL(I) ,FYR(I) 

Rotor elastic Y-force to the immediate left 
and right of a rotor station I, respectively 

pounds 

FYY 

A Y-force function 


GAL 

Combined shear rigidity 


IBI , INS 

Subscripts used in , dimensioned variables 


IR1 

A rotor section torsional rigidity variable 


KNS 

Subscript used in a dimensioned variable 


MDL,MDR 

Absolute value of the rotor transverse moment 
whirl to rotor spin velocity ratio at the 
immediate left and right of a rotor station, 
respectively 


MHX(I) ,MHY (I) 

Combined rotor transverse vis cons and Coulomb 
friction hysteresis XZ-plane and YZ-plane 
moment at rotor station I, respectively 

lb-in. 

MLS, MRS 

Rotor transverse elastic moment vector at the 
immediate left and right of a rotor station, 
respectively 

lb-in. 

MWL,MWR 

Rotor moment spin-whirl velocity difference to 
its absolute value ratio at the immediate left 
and right of a rotor station 

dimensionless 

MXL(I) ,MXR(I) 

Rotor XZ-plane elastic moment at the immediate 
left and right of rotor station I 

lb-in. 

MYL(I) ,MYR(I) 

Rotor YZ-plane elastic moment at the immediate 
left and right of rotor station I 

lb- in. 
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TABLE XX - (Continued) 


Variable 

Definition 

Units 

NB4 

Four times total number of bearings 


PPL 



SXL 

=SVXL+SCXL 


SXR 

=XVXR+SCXR 


SYL 

=SVYL+SCYL 


SYR 

=SVYR+SCYR 


THH 

T to the power HA 


T0L 

Tolerance used in limiting the computation 
round off error from generating unrealistic 
Coulomb friction induced hysteresis effect 

dimensionless 

T0R 

Net externally applied rotor drive and damping 
torque in an assumed rigid rotor sections 

lb-in. 

XBM(K) 

Bearing XZ-plane slope at bearing station K 

radians 

XMX 

Rotor X-displacement differential between 
the end stations of a rotor section 

inches 

XPL 



YBM 

Bearing YZ-plane slope at bearing station K 

radians 

YMY 

Rotor Y-displacement differential between the 
end stations of a rotor section 

inches 

YPL 



BCXL, BCXR 

Rotor Conlomb friction induced hysteresis 
bending moment in XZ-plane at the immediate 
left and right of a rotor station 

lb- in. 

BCYL, BCYR 

Rotor Conlomb friction induced hysteresis 
bending moment in YZ-plane at the immediate 
left and right of a rotor station 

lb-in 

1 
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TABLE XX - . (Continued) 


Variable 

Definition 

Units 

BVXL , BVXR 

Rotor viscous hysteresis coefficient induced 
XZ-plane bending moment at the immediate left 
and right of a rotor station 

lb-in. 

BVYL.BVYR 

Rotor viscous hysteresis coefficient induced 
YZ-plane bending moment at the immediate left 
and right of a rotor station 


COMB 

Combined torque 


EI1L,EI2L, 

EI3L 

1 

Bending sectional rigidity variables 


FD0D 

Rotor acceleration variable 


FWIL,FWIR 

The magnitude of rotor force spin-whirl 
velocity difference at the immediate left and 
right of a rotor station 

radians/sec 

FXLD 

Rotor force left X-velocity 

lb/sec 

FXRD 

Rotor force right X-velocity 

lb/sec 

FYLD 

Rotor force left Y-velocity 

lb/sec 

FYRD 

Rotor force right Y-velocity 

lb/sec 

IBNS 



I10S,I2NS, ) 
I3NS ,I4NS 
I5NS,I6NS 
I7NS,I8NS 
I9NS,K2NS 
K3NS , K4NS 
K5NS,K6NS 
K7NS J 


Subscripts 


MBVC 

MSVC 

MTVC 


Indicators for hysterisis input conditions 











TABLE XX - (Continued) 


Variable 

Definition 

Units 

MTXZ(I) ,MTYZ (I) 

Transverse torque XZ-plane and YZ-plane 
moment, respectively 

lb-in. 

MWIL,MWIR 

The magnitude of rotor whirl spin-whirl 
velocity difference at the immediate left 
and right of a rotor station 


MXLD(I) , 
"MXRD(I) 

Rotor elastic XZ-plane moment at the immediate 
left and right of the rotor station I 


MYLD(I) , 
MYRD(2) 

Rotor elastic YZ-plane moment at the immediate 
left and right of the rotor station I 


NST1 

An index in torsionally rigid rotor section 
control 


SCXL, SCXR 

Rotor Coulomb friction induced hysteresis shear 
force along X-axis at the immediate left and 
right of a rotor station 

pounds 

SCYL, SCYR 

Rotor Coulomb friction induced hysteresis shear 
force along Y-axis at the immediate left and 
right of a rotor station 

pounds 

i 

SQRT 

Square root library function 


SVXL , SVXR 

Rotor viscons coefficient induced hysteresis 
shear force along X-axis at the immediate left 
and right of a rotor station 

pounds 

SVYL.SVYR 

Rotor viscons coefficient induced hysteresis 
shear force along Y-axis at the immediate left 
and right of a rotor station 

pounds 

TMTX(I) , 
TMTY(I) 

Transverse torsion moment loading in XZ- and 
YZ-plane, respectively 

pounds 

T0RQ (J) 

Internal torque transmission including elastic 
torque hysteresis torque at rotor section J for 
computing the transverse loading effects of 
torsion 

lb-in. 

T0RS(I) 

Externally applied net rotor drive and damping 
torque at rotor station I 

lb-in. 

1 
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TABLE XX - (Continued) 


Variable 

Definition 

Units 

WHIR 

Rotor displacement whirl velocity 

radians/sec 

XF0R 

Rotor Y-force 

pounds 

XM0M 

Rotor XZ moment 

lb-in. 

YF0R 

Rotor Y-force 

pound 

YM0M 

Rotor YZ -moment 

lb-in. 

C0SFA 

Cosine of eccentricity angle 


C0SFG 

Cosine of misalignment angle 


EIC0M 

Rotor section rigidity 

lb-in. * 

FD0AB,FD0CT 

Spin velocity variable 

dimensionless , 
radians/sec 

FD0MT 

A spin-speed variable 

radians/sec^^ 

FD0RA 

A spin-speed ratio 

dimensionless 

FD0TM.FD0TN 

FD0TP 

Spin-speed differences 

radians/sec 

FXLD1,FXLD2 

Rotor X-force velocities at the immediate 

lb/sec 

FXRD1,FXRD2 

left and right of a rotor station, respectively 


FYLD1,FYLD2 

Rotor Y-force velocities at the immediate 

lb/sec 

FYRD1.FYRD2 

left and right of a rotor station, respectively 


IB2NS,IB3NS, 



IB4NS.IB5NS 
IB6NS, IB7NS 

Subscripts 


ISTAR 

Rotor torsional rigidity variables 


IST0P 



I10SB 

Subscript 


MXLD1,MXLD2 

Rotor XZ-plane moment velocities at the 

lb-in. /sec 

MXRD1,MXRD2 

immediate left and right of a rotor station 










TABLE XX - (Concluded) 


Variable 

Definition 

Units 

MYLD1 

MYLD2 

MYRD1 

MYRD2 

Rotor YZ-plane moment velocities at 
immediate left and right of a rotor 

the 

station 

lb-in. /sec 

NST0T 

Torsion rigidity variable 



SINFA 

Sine of eccentricity angle 



SINFG 

Sine of misalignment angle 



WHIRM 

Rotor moment whirl velocity 


^radians/sec 

XBD0T(I) 

Bearing X-velocity 


in. /sec 

XBF0R(I) 

Bearing X-force 


pounds 

XBMD0 

Bearing XZ-plane slope velocity 


radians/sec 

XBM0M 

Bearing XZ-plane moment 


lb-in. 

YBD0T 

Bearing Y-velocity 


in. /sec 

YBF0R 

Bearing Y-force 


pounds 

YBMD0 

Bearing YZ-plane slope velocity 


radians/sec 

YBM0M 

Bearing YZ-plane moment 


lb-in. 

FD0MAB 

FD0MRA 

FD0NAB 

FD0NRA 

Rotor spin velocity ratios 



FD0TSQ 

Square of rotor spin velocity 


radians/sec^ 

GALE I 3 

Rotor section rigidity variable 


in. /lb 

I10S2B,I10S3B 
I10S4B.I10S5B 
I 10S6B, I 10S7B 

Subscripts used in dimensioned variables 


TORHFM 

Hysteresis torque 


lb- in. 







• ^ APPENDIX C 

PROGRAM INPUT VARIABLES 

Appendix C consists of two sets of program input variables 

(a) NAMELIST/MUST - These variables must be input (Table XXI) 

(b) NAMELIST/0PTI0N - The input of these variables is optional, 
default values of the variables are as shown in Table XXII. 



TABLE XXI. PROGRAM INPUT VARIABLES (NAMELIST/MUST) 


Variables 

Description 

Unit 

Symbol Used 
and Eq, No* 

Title 

A description card consisting of 72 characters 
and ID field 




DT 

Initial integration step 

sec 




For Runge-Kutta and Adams-Moulton fixed-step 
techniques, DT will be used for the entire 
computation 





For predictor and corrector Adams-Moulton 
variable step technique, the input DT will 
be used only when the integration toler- 
ance requirement is met 




TMAX 

Maximum run real time 

sec 



DP 

Output printing real time interval 

sec 



NS 

Number of rotor stations of the rotor model 


"n" 

6-13 

NB 

Number of bearing stations 




FD0TI 

Initial rotor spin and whirl velocity 

r pm 

♦i 

1-5 

IB (I) 

Rotor station number for each bearing 




DD (J) 

Outside diameter for rotor sections 

cm 

D oi 

29-30 

QL(J) 

Rotor section length 

cm 

a. 

1 

20-29 

MET 

Unit system selector 




1 = international units 
0 = English units 




ID 

Not used in computation. It is a way to punch 
ID on cards for namelist input, so that the 
cards sequence can be defined 
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TABLE XXII - PROGRAM INPUT VARIABLES (NAMELIST/0PTI0N) 


Variable 

Description 

Default 

Value 

Unit 

Symbols Used 
and Eq. No. 

IND 

Integration technique selector 

0 = using Adams-Moulton predictor-corrector 

variable step technique 

1 = using 4th order Runge-Kutta fixed-step 

technique 

2 = using Adams-Moulton fixed-step 

technique 

1 



T0LI 

Computation accuracy control tolerance for Adams- 
Moulton variable step integration technique 

0.0001 



T 

Starting computation real time 

0 

sec 

t 28(i) ,46 

C0NTIN 

Specification for cold start or continued analysis 
from a previous un 

0 = cold start 

1 = continued analysis in which punch cards for 

T, DT, YNN(I) from a previous analysis must 
be included in addition to the input re- 
quired for cold start. T, DT value in the 
punched cards will be used to override 
that from the key punch input. 

0 



IT0RQ 

Input driving and damping torque control 

IT0RQ = 0 means no driving or damping torque 
effect will not be considered in 
computation 

IT0RQ = 1 means the torque effects will be 
considered in computation 
This is to bypass certain computation when 
IT0RQ = 0, to better computation efficiency 

0 



IPP 

Axial loading transverse effects control 
IPP = 1 including the effects 
IPP = 0 excluding the effects 

0 



IMT 

Torsional transverse effects control 
IMT = 1 including the effects 
IMT = 0 excluding the effects 

0 



RIG(J) 

Rotor section torsional flexibility indicator 
RIG(I) = 1 rotor section (I) is rigid in 
torsion 

RIG (I) = 0 rotor section (I) is torsionally 
flexible 

0 



CRT 

CRT graphs requirement specification 
CRT = 1 requires CRT 
CRT = 0 no CRT required 

0 



M0SHA 

Rotor mode shape CRT generation control 

1 



*NP0INT 

Number of points per CRT graph range from 
1 to 50 

25 



N00RPM 

Number of rotor spin speeds one for each rotor 
mode shape CRT plot 

1 



IASIGN 

Rotor station number for which the related CRT 
graph will be plotted 

1 



INPRPM(M) 

The values of input rotor spin speeds at or 
near which rotor mode shape CRT will be 
plotted 

0 

rpm 


D(J) 

Rotor section inside diameter 

0 

cm 

D u 29-30 
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TABLE XXII - (Continued) 


Variable 

Description 

Default 

Value 

Unit 

Symbols Used 
and Eq. No. 

DN(J) 

Rotor mass density 

0.008304 

kg/cm 3 


P(J) 

Rotor section Poisson’s ratio 

0.3 

' 


EE(J) 

Rotor section Young's modulus of elasticity 

2.0684 
x 10 7 

2 

Newtons/cm 

E. 20-23,36- 
1 39,42-45 

GG(J) 

Rotor section shear modulus of rigidity 

0.7929 
X 10 7 

Newtons/cm 

G 20,21,34, 
35,40,41 

EI(J) 

Product of rotor modulus of elasticity and area 
moment of inertia 

0 

Newton/cm 2 


GAK(J) 

Product of rotor shear modulus, cross-sectional 
area and shear strain over average shear 
strain ratio 

0 

Newtons 


AM(I) 

Additional (nonstructural) rotor masses 

0 

kg 


ECC(I) 

Rotor mass eccentricity 

0.000254 

cm 

e i 1,2,3 

ALFA (I) 

Rotor mass eccentricity phase angle 

0 

degrees 

a. 1,2,3 

AIR0 (I) 

Additional (nonstructural) rotor polar mass 
moments of inertia 

0 

2 

kg -cm 


AID (I) 

Additional (nonstructural) rotor transverse 
mass moments of inertia 

0 

2 

kg-cm 


BETA(I) 

Rotor inertia axis misalignment 

0 

degrees 

4,5 

GAMMA(I) 

Rotor inertia axis misalignment phase angle 

0 

degrees 

Y t 4,5 

BKMX(K) 

Mount in-phase anisotropic stiffness force 

3.5025 

N/cm 

^ 17a 
*wyi 

BKMY(K) 

coefficient along x- and y-axis, respectively 

xlO 6 

each 


BCMX(K) 

Mount in-phase anisotropic damping force coef- 

0 

(N-cm-sec)/ 

r Mxl 17b 
Myi 

BCMY(K) 

ficients along x- and y-axis, respectively 


cm 

XKMM(K) 

Mount in-phase anisotropic stiffness moment 

22.59697 

(N-an)/ 

K AMvi 

XKMM(K) 

coefficients in xz- and yz-plane, respectively 

x 10 6 
each 

radian 

* Myl 17c 

Vlyi 

XCMM(K) 

Mount in-phase anisotropic damping moment coef- 

22.59697 

(N-cm-sec)/ 

^ 17d • 

4>Myi 

YCMM(K) 

ficients in xz- and yz-plane, respectively 

x 10 6 
each 

radian 

BM(K) 

Bearing mass 

0 

kg 

M b . 14,15 

BI(K) 

Bearing transverse mass moment of inertia 

0 

2 

kg-cm 

I Bi 16.17 

QKXX(K) 

Bearing in-phase anisotropic stiffness force 

1.7513 

Newtons/cm 

H*** 1 3, 
Byyi 

QKYY(K) 

coefficients along x- and y-axis, respectively 

x 10 6 
each 


QKXY(K) 

Bearing out-of-phase anisotropic stiffness 
x-force due to y-displacement coefficient 

0 

Newton/ cm 

K d . 3a 
Bxyi 

QKYX(K) 

Bearing out-of -phase anisotropic stiffness 
y-axis due to x-displacement coefficient 

0 

Newton/cm 

Vi 3a 

QCXX(K) 

Bearing in-phase anisotropic damping force 

0 

(Newton-sec)/ 

r Bxxi 3a 
C Byyi 

QCYY(K) 

coefficients along x- and y-axis, respectively 

I 

cm 
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TABLE XXII - (Continued) 


Variable 

Description 

Default 

Value 

Unit 

Symbols Used 
and Eq. No. 

QCXY(K) 

Bearing out-of-phase anisotropic damping 
x-force due to y~velocity coefficient 

0 

(Newton-sec)/ 

cm 

C Bxyi 3a 

QCYX(K) 

Bearing out-of-phase anisotropic damping 
y-force due to x-velocity coefficient 

0 

(Newton-sec)/ 

cm 

C Byxi 3a 

XXMK(K) 

YYMK(K) 

Bearing in-phase anisotropic stiffness moment 
coefficients in xz- and yz-plane, respectively 

11.29848 

x 10 6 
each 

(Newton-cm)/ 

radian 

> xxi 5a 
B$yyi 

XYMK(K) 

Bearing out-of-phase anisotropic stiffness 
xz-plane moment due to yz-plane slope 
coefficient 

0 

(Newton-cm)/ 

radian 

K B<J>xyi 5b 

YXMK(K) 

Bearing out-of-phase anisotropic stiffness 
yz-plane moment due to xz-plane slope 
coefficient 

0 

(Newton-cm)/ 

radian 

Vxi 5b 

XXMC(K) 

YYMC(K) 

Bearing in-phase anisotropic damping moment 
coefficients in xz- and yz-plane, respectively 

0 

(Newt on- cm - 
sec) /radian 

%«xi Sc 
B4>yyi 

XYMC(K) 

Bearing out-of-phase anisotropic damping 
xz-plane moment due to yz-plane slope velocity 
coefficient 

0 

(Newton-cm- 
sec) /radian 

C n . . 5d 
B<J>xyi 

YXMC(K) 

Bearing out-of-phase anisotropic damping 
yz-plane moment due to xz-plane slope velocity 
coefficient 

0 

(Newton-cm-. 
sec) /radian 

C D . . 5d 
B<{>yxi 

KK(K) 

Number of nonlinear bearing stiffness sections 

1 



FD0FIX(I) 

Nonlinear stiffness bearing rotor spin-speed 
factor 

0 

radians/sec 

*oi 2 * 3 

BBB(K,K1) 

Nonlinear stiffness bearing spin-speed, 
journal displacement coefficient 

0 

(Newton-sec)/ 

radian 

B Bik 2 ' 3 

BCB(K.Kl) 

I 

Nonlinear stiffness bearing journal displace- 
ment power coefficient 

0 

l/(cra) BHB (K. K1 ) 

C Bik 2 ' 3 

BDB(K,K1) 

Nonlinear stiffness bearing journal displace- 
ment coefficient 

0 

cnf * 

°Bik 2 ’ 3 

BEB(K.Kl) 

Nonlinear stiffness bearing constant 

0 

dimensionless 

E Blk 2 ' 3 

BHB(K,K1) 

Nonlinear stiffness bearing displacement 
exponent 

1.0 

dimensionless 

H Bik 2 « 3 

BKB(K,K1) 

Nonlinear stiffness bearing coefficient 

0 

Newtons/cm 

K Bik 2 ' 3 

BNB(K,K1) 

Nonlinear stiffness bearing spin-speed 
coefficient 

0 

(Newton- sec)/ 
(cm-radian) 

N Blk 2 > 3 

BR0B(K,KZ) 

Nonlinear stiffness bearing stiffness lower 
limit bearing displacement for a stiffness 
section 

BR0B (K, 1) 
=0 

BR0B (K,2) 
= .0127 

cm 

P 2 3 

Bik 

QK(I) 

Rotor-to-casing in-phase stiffness force 
coefficient 

0 

Newtons /cm 

K. 2,3 

QC Cl) 

Rotor-to-casing in-phase damping force 
coefficient 

0 

(Newton-sec)/ 

cm 

C. 2,3 

QKP(I) 

Rotor-to-casing out-of-phase stiffness force 
coefficient 

0 

Newtons/cm 

K . 2,3 
Pi 
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TABLE XXII- (Continued) 


Variable 

Description 

Default 

Value 

Unit 

Symbols Used 
and Eq. No. 

QCP(I) 

Rotor-to-casing out-of-phase damping force 
coefficient 

0 

(Newton-sec)/ 

cm 

C . 
Pi 

2,3 

QKF(I) 

Rotor-to-casing in-phase stiffness moment 
coefficient 

0 

(Newton-cm)/ 

radian 

K.. 

4,5 

QCF(I) 

Rotor-to-casing in-phase damping moment 
coefficient 

0 

(Newton-cm- 
sec) /radian 

C . . 
4>i 

4,5 

QKPF(I) 

Rotor- to-casing out-of-phase stiffness moment 
coefficient 

0 

(Newton-cm)/ 

radian 

K. . 
(frpi 

4,5 

QCPF(I) 

Rotor- to- casing out-of-phase damping moment 
coefficient 

0 

(Newton-cm- 
sec) /radian 

C< tpi 

4.5 

XKF(I) 

Rotor-to-casing whirl stiffness force factor 

0 

dimensionless 

K Fi 

2,3 

XCF(I) 

Rotor-to-casing whirl damping force factor 

0 

dimensionless 

C Fi 

2,3 

XKFF(I) 

Rotor-to-casing whirl stiffness moment factor 

0 

dimensionless 

K. u . 

4>Mi 

4.5 

XCFF(I) 

Rotor- to-casing whirl damping moment factor 

0 

dimensionless 

C <t>Mi 

4.5 

QKHD(I) 

Rotor-to-casing out-of-phase whirl-spin stiff- 
ness force coefficient 

0 

(Newton-sec)/ 

cm 

*SlDi 

2,3 

QCHD(I) 

Rotor-to-casing out-of-phase, whirl-spin damping 
force coefficient 

0 

(Newton-sec 2 )/ 

radian 

C HDi 

2,3 

QKHDF(I) 

Rotor-to-casing out-of-phase, whirl-spin stiff- 
ness moment coefficient 

0 

(Newton-cm- 
sec) /radian 

K <f)HDi 

4,5 

QCHDF(I) 

Rotor- to-casing out-of-phase, whirl-spin damping 
moment coefficient 

0 

(Newt on- cm- 
sec 2 ) /radian 

C <t>HDi 

4,5 

CT1(I) 

Rotor damping torque spin- speed power coefficient 

0 

(Newton-cm) (sec/ 
radian) 01 ^ 

C Tli 

28 

CT(I) 

Rotor damping torque spin- speed 

0 

dimensionless 

C Ti 

28 

CI2 (I) 

Rotor damping torque spin-speed coefficient 

0 

(Newton- cm- 
sec) /radian 

C T2i 

28 

HT(I) 

Rotor drive torque spin-speed exponent 

0 

dimensionless 

*Ti 

28 

MT1(I) 

Rotor drive torque spin-speed power coefficient 

0 

(Newton-cm) (sec/ 
radian)^ 1 ) 

"ru 

28 

MT2 (I) 

Rotor drive torque spin-speed coefficient 

0 

(Newt on- cm - 
sec) /radian 

k 2i 

28 

AT(I) 

Constant rotor drive torque 

0 

Newton-cm 

^Ti 

28 

BT(I) 

Rotor drive torque time coefficient 

0 

(Newton-cm)/ 

sec 

B n 

28 

DU(I) 

Rotor drive torque time power coefficient 

0 

(Newton-cm)/ 

(sec) 1 ^ 1 ) 

D Ui 

28 

ET(I) 

Rotor drive torque sine coefficient 

0 

Newton-cm 

E Ti 

28 

HT(I) 

Rotor drive torque time power exponent 

0 

dimensionless 

«ri 

28 

FT(I) 

Rotor drive torque time coefficient sine function 

0 

radians/sec 

F Xi 

28 
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TABLE XXII - (Concluded) 




Default 


Symbols Used | 

Variable 

Description 

Value 

Unit 

and Eq 

. No. 

GT(I) 

Rotor drive torque constant argument for sine 
function 

0 

radians 

G Ti 

28 

AA(I) 

Rotor axial loading constant 

0 

Newtons 

\i 

46 

BA(I) 

Rotor axial loading time coefficient 

0 

Newtons/sec 

B Ai 

46 

DA(I) 

Rotor axial loading time power coefficient 

0 

Newton/ (sec) ^ 

°Ai 

46 

EA(I) 

Rotor axial loading sine coefficient 

0 

Newtons 

e a . 

Ai 

46 

FA 

Rotor axial loading sine function time coefficient 

0 

radians/sec 

FA 

46 

HA 

Rotor axial loading time power exponent 

0 

dimensionless 

HA 

46 

GA 

Rotor axial loading since function constant 
argument 

0 

radians 

GA 

46 

GX 

Acceleration or gravity loading along negative 

0 

2 

cm/sec 

g x 

2 

GY 

x- and y-axis, respectively 

0 



3 

USV(J) 

Rotor transverse shear viscous hysteresis 

0 

(Newton- sec) / 

P SVi 

34,35 


coefficient 


2 

cm 


USC(J) 

Rotor transverse shear Coulomb friction hysteresis 
coefficient 

0 

Newtons/cin 

y SCi 

40,41 

UBV(J) 

Rotor transverse bending viscous hysteresis 

0 

(Newton- sec) / 

y BVi 

36-39 


coefficient 


cm 2 


UBC(J) 

Rotor transverse bending Coulomb friction 
hysteresis coefficient 

0 

2 

Newtons/cm 

p BCi 

40-43 

UTV(J) 

Rotor torsional shear viscous hysteresis 

0 

Newton-sec/ 

u TVi 

29 


coefficient 


2 

cm 


UTC(J) 

Rotor torsional shear Coulomb friction coefficient 

0 

2 

Newtons /cm 

y TCi 

30 

FI 

Initial rotor spin angular position 

1.0 x 

degrees 

♦i 

2-5 



10- 20 



IPRINT 

IPRINT x DP = computer output printout real 
time interval 

1 

sec 



ID 

Not used in computation, it is a way to punch ID 
on cards for namelist input, so that the card 
sequence can be defined 
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APPENDIX D 


A a . 

Ai 

A'. 

1 

b. . 
il 


B Ai 

B BiK 


B Ti 

C i*Si 


C. . 


'BiK 


Bxxi 

C D . 

Byyi 

c D . 

Bxyi 

C n ' . 

Byxi 

^B<J>xxi 

^B(j)yyi 

^B(}>xyi 

^Bijiyxi 

Slxi ’Slyi 

C Pi ,C <J>Pi 


C Ti 


"Tli 


"T2i 


SYMBOLS LIST FOR THE MATHEMATICAL FORMULATION 

rotor axial loading constant force. Newtons 

rotor cross-sectional area, cm 

rotor constant drive torque, Newton-cm 

rotor displacement (j ) -moment (i) influence coefficient, 
cm/ (Newton-cm) 

rotor axial loading time coefficient, Newtons/sec 

nonlinear stiffness bearirtg spin-speed and journal displacement 
coefficient. Newton-sec/radian 

rotor drive torque time coefficient, (Newton-cm) /sec 

in-phase force and moment damping coefficient, respectively, 
(Newton-sec) / cm, (Newton-cm-sec) /radian 

rotor displacement (j) -force (i) influence coefficient, 
cm/Newton 

nonlinear stiffness bearing displacement power coefficient, 
(l/cm)" 81 * 

in-phase anisotropic bearing damping force coefficients along 
x- and y-axis, respectively, (Newton-sec) /cm 


out-of-phase anisotropic bearing damping coefficients along 
x- and y-axis due to velocities along y- and x-axis, respectively 
(Newton-sec) /cm 

in-phase anisotropic bearing damping moment coefficients in 
xz- and yz-plane, respectively, (Newton-cm-sec) /radian 

out-of-phase anisotropic bearing damping moments in xz- and yz- 
plane due to rotations in yz- and xz-plane, respectively, 
(Newton-cm-sec) /radian 

mount in-phase damping force coefficients along x- and y-axis, 
respectively, (Newton- sec). /cm 

out-of-phase force and moment damping coefficients, respectively, 
(Newton-sec) /cm, (Newton-cm-sec) /radian 

rotor torsion spin-speed damping spin-speed power coefficient 
exponent, dimensionless 

rotor torsion spin-speed damping spin-speed power coefficient, 
(Newton-cm) (sec/radian) 

rotor torsion spin-speed damping coefficient, (Newton-cm-sec)/ 
radian 
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"HDi 

: 4>hd. 


C Fi ,C (J>Fi 


4>Mxi 

G 4>Myi 

°Ai 

°BiK 

D Ii ,D oi 

D Ui 


e. 

i 


Ai 

E. 

1 

E BiK 

E Ti 


Ti 


F . ,F . 
xi* yi 

g X »Sy 


G. 

1 

G Ti 


H A 

H BiK 

HL,HR 


HSC,HSV 

HBC,HBV 

“li 


out-of-phase whirl and spin speed sensitive force and moment 

2 

damping coefficients, respectively, (Newton-sec ) / (cm-radian) , 

2 2 

(Newton-cm-sec ) /radian 

rotor-to-casing whirl damping force and moment factor, respectively, 
dimensionless 

mount in-phase damping moment coefficients in xz- and yz-plane, 
respectively, (Newton-dm-sec) /radian 

rotor axial loading time power coefficient, Newton/ (sec)^ 
nonlinear stiffness bearing journal displacement coefficient, 1/cm 
rotor section inside and outside diameter, respectively, cm 

Hri 

rotor drive torque time power coefficient, (Newton-cm)/(sec) 11 
rotor mass eccentricity, cm 

rotor axial loading sinussoidal function coefficient. Newtons 

2 

rotor Young's modulus of elasticity. Newtons/cm 

nonlinear stiffness bearing constant, dimensionless 

rotor drive torque sinusoidal function coefficient, Newton-cm 

rotor axial loading sinusoidal function constant argument, radians 

rotor drive torque sinusoidal function argument time coefficient, 
radians/sec 

force loading on rotor along +x- and +y-axis, respectively, Newtons 
gravitational or g-loading along -x- and -/-axis, respectively, 
cm/sec^ 

rotor shear modulus of rigidity, Newton/cm^ 

rotor drive torque sinusoidal function constant argument, radians 

axial loading sinusoidal function constant argument, radians 

rotor axial loading time power exponent, dimensionless 

nonlinear stiffness bearing displacement exponent, dimensionless 

subscripts pertaining to rotor hysteresis force or moment at 
left or right of a rotor station 

subscripts pertaining to the shear Coulomb friction and viscous 
hysteresis force coefficient, respectively 

subscripts pertaining to the bending Coulomb friction and viscous 
hysteresis moment coefficient, respectively 

rotor drive torque time power exponent, dimensionless 
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1 

I. 


Bi 

I Di ,I Pi 


Bxxi 

C Byyi 

Sxyi 

Syxi 

S<f>xxi 

Scfiyyi 

S<f>xyi 

Stfiyxi 


K BiK 

K Fi 

W 


HDi 

K 4»HDi 

K Ii 

K,. r 

<J)Mi 

^Mxi 

*Slyi 

K Pi 

K $Pi 

% 

^4>Mxi 

*Sj>Myi 


subscript pertaining to ith rotor station or section as appropriate 

2 

rotor area moment of inertia, cm 
bearing transverse mass moment of inertia, kg-cm^ 
equivalent discrete rotor transverse and polar mass moment of 
inertia, respectively, kg-cm^ 

subscript pertaining to kth stiffness section for nonlinear 
stiffness bearing 

in-phase anisotropic bearing stiffness force coefficients along 
x- and y-axis, respectively. Newton/cm 


out-of-phase anisotropic bearing stiffness force along x- and y- 
axis due to displacements coefficients along y- and x-axis, 
respectively, Newtons./cm . 

in-phase anisotropic bearing stiffness moment coefficients in 
xz- and yz-plane, respectively, (Newton-cm) /radian 

out-of-phase anisotropic bearing stiffness moment in xz- and yz- 
phase due to yz- and xz-plane rotations coefficients, respec- 
tively, (Newtons-cm) /radian 

nonlinear stiffness bearing coefficient, Newton/cm 

rotor-to-casing whirl stiffness force factor, dimensionless 

in-phase force and moment stiffness coefficient, respectively, 
Newton/cm, (Newton-cm) /radian 

out-of-phase force and moment stiffness coefficient, respectively, 
(Newton-sec) / (cm-radian) , (Newton-cm-sec) /radian^ 


rotor-to-casing whirl stiffness force and moment factor, respec- 
tively, dimensionless 


mount in-phase stiffness force coefficients along x- and y-axis, 
respectively, Newtons/cm 


out-of-phase force and moment stiffness coefficients, respectively, 
Newtons/cm 

rotor sectional torsional stiffness, (Newton-cm) /radian 

mount in-phase stiffness moment coefficients in xz- and yz-plane, 
respectively, (Newton-cm) /radian 
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1. 

1 

L 

M. 

4 

M . 

XZl 


M 


yzi 

/ 


M 


xzi 


M 


yzi 

*ri 

^rii 




N BiK 


Q 

s 

s', .s'. 

xi ' yi 


Fij 


Mi j 


X. , Y . ,Z, 

ill 


Wi 

VV z i 

x bi ,y bi 


X Bi ,Y Bi 


X Bi ,Y Bi 

*bNB 

Y 

dNB 


= rotor section length between stations i and i+1, cm 
= z-distance between first and last bearing stations, cm 
= equivalent discrete rotor mass, kg 
= bearing mass, kg 

= moment rotor loading in x.z. and y.z. plane, respectively, 
Newton-cm 11 11 

= rotor transverse elastic moment as, defined in Fig. 5 


= rotor torsional drive torque spin-speed power exponent, 
dimensionless 


= rotor torsional drive torque spin-speed power coefficient, 
Newton-cm- (sec/radian) Ti 

= rotor torsional drive torque spin-speed coefficient, (Newton-cm- 
sec) /radian 

= nonlinear stiffness bearing spin-speed coefficient, (Newton-sec)/ 
(cm-radian) 

= z-coordinate of the last bearing station, cm 

= z-coordinate of the first bearing station, cm 

= rotor transverse elastic shear as defined in Fig. 5 

= rotor slope rotation (j) -force (i) influence coefficient, radian/ 
Newton 


= rotor slope rotation (j) -moment (i) influence coefficient, radian/ 
(Newton-cm) 


= a right-hand rotor system displacement coordinates with and 
defined in Fig. 

= first time derivations of X., Y. , Z. , cm/sec 

ill 2 

= second time derivatives of X^, Y^, Z^, cm/sec 

= bearing displacements at the first bearing station along x- 
and y-axis, respectively, cm 


= journal x- and y- displacement with respect the bearing, respec- 
tively, cm 


= first time derivatives of X n ., 


Yg^, respectively, cm/sec 


= bearing displacement at the last bearing station along x- and y- 
axis, respectively, cm 
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*Mi ’ Y Mi 

MmI " 

V V Mi = 

Z. 

l 

Greek Symbols 


mount displacements along x- and y-axis, respectively 
first time derivatives of x^ and Y^, respectively, cm/sec 
second time derivatives of and Y^, respectively, cm/sec 
z-coordinate of rotor station i, cm 


2 


a. 

1 


a. 

l 

6 i 


Bxzi 


Byzi 

*Bxzi * 
*Byzi 


yzi 


.yzi 

3 xzi’ 

3 . 
yzi 

J BCi 

J BVi 


SCi 1 

J SVi 


y TCi’ 

y TVi 

P BiK 

Mi 


= phase angle for e^, radians 

= shear deflection correction factor, dimensionless 

= rotor mass moments of inertia misalignment, radians 

= rotor mass moments of inertia misalignment phase angle, radians 

= journal slopes with respect to the bearing in xz- and yz-planes, 
respectively, radians 

= first derivatives of 0 D . and 0 D . , respectively, radians/sec 

bxzi byzi 


= angular displacements of rotor elastic centerline in xz- and yz- 
plane, respectively, radians 

= first time derivatives of 0 . , 0 . , radians/sec 

xzi yzi’ 


2 

= second time derivatives of 0 . , 0 . , radians/sec 

xzi’ yzi’ 

= rotor transverse bending Coulomb friction and viscous hysteresis 
coefficient, respectively, (Newton-sec) /cm2 , Newtons/cm2 


= rotor transverse shear Coulomb friction and viscous hysteresis 
coefficient, respectively, (Newton-sec)/cm2 , Newtons/cm2 


= rotor torsional shear Coulomb friction and viscous hysteresis 

coefficient, respectively, Newtons/cm^, (Newton-sec)/ (radians-cm2) 

= nonlinear stiffness bearing initial displacement, cm 

= rotor spin angular displacement, velocity, and acceleration, 
respectively, radians, radians/sec, and radians/sec^ 



nonlinear bearing speed sensitive coefficient, radians/sec 


W Fi ,t0 Mi = force and moment whirl velocity, respectively, radians/sec 
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APPENDIX E 


COMPUTER PROGRAM LISTING 

The rotor dynamics computer program listing, including the CDC system graphic 
plotting portions, is herewith attached. The graphic plotting portions in MAIN 
and in HYSMET, have not been checked out in accordance with contractual agree- 
ments. 

Five sets of namelist data used in obtaining the final computer verification 
are also attached for reference. 
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APPENDIX F 


IBM 360/370 COMPUTER RESULTS 

Tabled 111 provides the computation results. for a transient spin speed 
rotor operation in the International Unit System. Table XXIV is the con- 
stant speed rotor dynamic performance for a rotor configuration having 15 
rotor stations and 6 support bearings, which represents the current maxi- 
mum computer program capacity. 
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TORSIONAL SHEAR VISCOUS COEFFICIENT ARRAY (UTV(J)), LE*SFC/IN**2 TABLE XXIV. Continued 
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